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ABSTPvACT 

We develop and implement an observational test of the theoretical notion that dissipation in major mergers 
of gas-rich galaxies produces the fundamental plane (FP) and related correlations obeyed by ellipticals. Obser- 
vations have shown that the "tilt" of the FP involves more than a simple non-homology or stellar population 
effect: lower-mass ellipticals have a higher ratio of stellar to dark matter within their stellar effective radii. 
Theoretical models have attempted to explain this via dissipation: if ellipticals are formed in major mergers 
of disks, then mergers between disks having a larger gas content (typically observed to be lower-mass disks) 
will yield remnants with a larger mass fraction formed in a central, compact starburst, giving a smaller stellar 
R e and lower M tot /M* within that R e . Such starbursts leave a characteristic imprint in the surface brightness 
profiles of ellipticals, in the form of a central excess above the outer profile established by the dissipationless, 
violent relaxation of disk stars. In previous work, we implemented a purely empirical method to use such 
features in the observed profiles of ellipticals to robustly estimate the amount of dissipation involved in the 
original spheroid-forming merger. Applying this to a large sample of ellipticals with detailed kinematic and 
photometric observations, we demonstrate that the location of ellipticals on the FP and its tilt are in fact driven 
by dissipation. We show that at fixed mass, ellipticals formed in more dissipational events, as indicated by their 
observed profiles, are smaller and have a lower ratio M tot /M*. At the same (fixed) degree of dissipation, there 
is no tilt in the FP - i.e. ellipticals formed with a similar level of dissipation have the same ratio of enclosed 
stellar to total mass within R e . 

We further demonstrate that observations and these models obey the "homology assumption," i.e. that the 
true enclosed mass M tme (R e ) oca 2 R e . Measured at the radii of disks of the same mass, we show that ellipticals 
have the same total enclosed masses as those disks - i.e. that the FP tilt can be effectively removed. Therefore, 
the fundamental plane tilt cannot primarily owe to non-homology or to changes in the dark matter distribution: 
it must arise as a result of a contraction of the baryonic component relative to the dark matter in the process that 
transforms disks to ellipticals, as predicted by dissipational mergers. If we allow for the observed cosmological 
dependence of disk gas fraction on mass, the observed FP, size-mass, and velocity dispersion-mass correlations 
are reproduced by our models, as are the observed homology constraints and profile shapes. Dissipation is both 
necessary and sufficient to explain the observed FP correlations of ellipticals. These observations all favor 
theories in which ellipticals are formed in major mergers of disks with gas fractions, sizes, and dark matter 
content similar to that observed as a function of mass in low-redshift disks: unusually compact disks are not 
required to make ~ 0.01 - 10L* ellipticals. We present a number of associated predictions that can be used to 
further test these assertions. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation — 
galaxies: nuclei — galaxies: structure — cosmology: theory 
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1. INTRODUCTION 

Understanding the scaling relations between the photo- 
metric and kinematic properties of galaxy spheroids - their 
masses, sizes, velocity dispersions, and luminosities - is 
fundamental to explain ing the origin of early-type galaxies. 
iFaber & Jackson! d 19761) demonstrated that ellipticals obey a 
relatively tight correlation bet ween optica l lumi nosity and 
central velocity dispersion, and iKormendyl (1 19771) found an 
analogous relationship between their effective surface bright- 
ness and radii. With improved observations and the advent 
of stellar population modeling, these observed trends can be 
translated into robust correlations between physical parame- 
ters: a velocity dispersion-stellar mass (er— M*) and a size- 
st ellar mass (R e -M*) r e lation. 

IDjorgovski & Davis! (119871) and iDressler et all (1 19871) 
demonstrated that the scatter in either the Fabe r & Jackson! 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138 

2 W. M. Keck Postdoctoral Fellow at the Harvard-Smithsonian Center for 
Astrophysics 



d 19761) or [Kormendy] d 19771) relation could be reduced by 
adopting a three parameter correlation of the form log(R e ) = 
a log (<t)- 0.4 b/j, e + c (equivalently R e cx <7 a /*), with best fit 
scalings a ^ 1.3—1.4, & ~ -0.8 to -0.9. This defines the "fun- 
damental plane" (FP) of elliptical galaxies: a correlation relat- 
ing stellar mass or luminosity (implicit in the surface bright- 
ness), effective radius, and velocity dispersion (effectively the 
dynamical mass of the system). With a small observed scatter 
~ 0.1 dex, the FP has presented as a long-standing, and still 
unresolved challenge to observations and theoretical models 
of spheroid formation. 

In developing a physical under standing of the FP an d 
associated ellip t ical scaling laws, IDjorgovski et al.l (Il988h : 
iJorgensen et al.l (1 19931) and others demonstrated that the FP 
could be represented as a "tilted" virial plane. If ellipti- 
cals were perfectly homologous systems with constant stel- 
lar mass-to-light ratios M*/L, then a virial correlation L oc 
= ka 2 R e /G = Mdyn, with constant integral factor k, would 
be expected. Since I e cx L/R 2 e , this translates to an expected 
"virial FP" R e cx a 2 1' 1 . The observed FP is similar to this, but 
not exactly so; it is equivalent to and can be represented as a 



2 



Hopkins et al. 



'tilted" version of this correlation, namely 



M iyn otMl +a (1) 

with some small, but non-zero a. Equivalently, the difference 
between the best-fit observed FP (in any projection) and the 
virial FP can be expressed as a mass dependent mass-to-light 
(or, for our purposes, total mass-to-stellar mass) ratio 
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where the quantity a quantifies the tilt, or deviation of the 
FP from the virial relation. Various indep endent measure- 
ments find similar values of a w 0.2 (e.g. [P ahre et al ] 1 19981; 
Gerhard et al. 2001; Borriello et al. 2003c IPadmanabhan et al.1 
2004tlGallazzi et al.ll2006h . Although this is not strictly iden- 
tical to the best-fit relation R e oc a" 7* if both a and b are fit 
as free parameters, multiple observations have shown that it 
is statistically an equivalent representation (i.e. has the same 
scatter in physical quantities), and that there is no additional 
information in the best-fit FP beyond this tilt (i.e. once this 
tilt is accounted for, there is no additional systematic scaling 
in R e or a that can reduce the scatter in predicting the other 
quantities). 

It is now well-established that part of the observed tilt in 
optical bands is a consequence of stellar population effects: 
lower-mass ellipticals tend to b e younger, yielding lower stel- 
lar mass-to-light ratios (see e.g. Trageret al. 2000). However, 
various constraints imply that only a small fraction of the 
optical tilt owes to these effects (see e.g. [P ahre et al. 1 1998 1 : 
[ Gerhard et all 120011; i Bertin et ail 120021: [Bordello et al J 120031 
Padmanabhan et al. 2004; Truiillo et al. 2004; Gallaz zi et al.l 
2006; von der Linden et aU 120071) . For example, in the K- 
band, the tilt is still substantial: observations indicate Md yn oc 
L L25±o.05 (Pa hre e t al] 1 19981) . whereas the systematic de- 
pendence of M*/Lk is quite weak (most estimates suggest 
M„/L K oc L^ 03 : lBelletal.ll2003h . It is now possible to com- 
bine high-resolution spectra and stellar population synthesis 
models, allowing reliable stellar mass estimates, and almost 
all such studies yield a similar relation: 
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i.e. a w 0.2 as described above. It has been demonstrated 
that this result is robust to the details of the stellar popula- 
tion model, spectral coverage, or even simplifying assump- 
tions such as the use of a single color to derive a mean M*/L, 
There are only two reasonable explanations for this find- 
ing (some combination of the two is also possible). First, 
the true mass enclosed within the stellar effective radius 
(M tme ) could in fact be proportional to the stellar mass M*, 
but owing to e.g. changes in the profile shape or kinemat- 
ics of galaxies with mass (traditional non-homology), the 
relation between actual mass and the dynamical mass esti- 
mator Md yn oc a 2 R e /G is a changing function of mass. In 
other words, Md yn oc M^ e . However, observations appear 
to rule out this possibility, at least as the origin of most of 
the tilt. Integral modeling of the mass distribution from two- 
dimensional kinematic maps (which should recover any sys- 
tematic difference between Md yn a nd without refere nce 
to any homology assumptions, e.g. Cappellari et al. 20061). as 
well as mass distributions estimated from gravitational lens- 



ing dBolton et all 12007. 2008: iNipoti etlLll2008l) . indepen- 
dently give M d yn oc M^ ±om . That is, the allowed contri- 
bution of non-homology to the FP tilt is small. 



The only remaining explanation is that the FP tilt reflects a 
meaningful physical change, namely that the ratio of total en- 
closed mass within R e (M tme , represented reasonably well to 
within a normalization constant by Md yn ) to the stellar mass is 
an increasing function of mass (M true oc M' 2 ). In other words, 
low-mass ellipticals are more baryon-dominated within their 
stellar R e , and high-mass ellipticals have higher dark matter 
fractions. We emphasize that these constraints apply within 
the stellar effective radii. The change in dark matter fraction 
is not required to be global: if one were to contract the stel- 
lar mass distribution but keep the dark matter halo relatively 
fixed, for example, it would significantly decrease the dark 
matter fraction (and correspondingly Mdyn/M») within R e . 

This trend is contrary to that followed by disks. For 
disk galaxies, an opposite (negative) t ilt Md vn oc m2' 7 ~ 0,8 
(a ~ -0.2 to -0-3) is observed (see e.g. iPersic et al.|[ l996b; 
iBell & de Jong|l200lt ISh^t^l200lTcourteau et al.H2007T) 
- low mass disks (and dwarf spheroidals) are the mo st dark- 
matter dominated systems (IPersic & Salucci| [19881 119901: 
IPersic et al.lll996albl:Bo"rrieUo & Saluccill200ll) . Such a scal- 
ing is expected if the properties of disks track those of 
their dark matter halos: lower-mass halos are more compact 
dNeto et al.ll2.007l and references therein), and it is also well- 
established that lower-mass disks experience less effici ent star 
formation (IBell & de JongH2000t iGallazzi et aHl2005h . Con- 
sequently, disks and ellipticals have similar ratios Mdy n /M» at 
high masses (~a few L» ), but disks are more dark-matter dom- 
inated (have much higher M& yn {R e )) at l° w (stellar) masses. 
This difference in scaling laws also relates to their stellar size- 
mass relations: disks obey a shallow relation R e oc M° - 25 " 035 , 
roughly consistent with the scaling of halo effective radii as 
a function of mass, whereas sphero ids obey a much steeper 
relation R e oc M° 6 (IShen et al.ll2003l) . Again, at -a few L*, 
disks and ellipticals have similar sizes and densities, but at low 
masses (<^_ L*), ellipticals are smaller (in their stellar/baryonic 
mass distributions) and more dense. 

This difference has, for ~ 30 years, represented a major 
challenge for theory - especially models which posit that 
ellipticals are formed thr ough the merger of disk galaxies 
(the " merger hypothesis"; iToomre & Toomrelll972t iToomrd 
1977). In particular, it has been noted that purely dissipation- 
less mergers of stellar disks cannot raise the mass and phase- 
space densities of ellipticals, and so cannot change th e scaling 
laws of ellipticals to something di fferent from disks (lOstrikerl 
[19801 ICarlberdfT986l lGurJl987h . 

However, these arguments do not pertain if ellipticals are 
formed from mergers of gas- rich disks. Particularly at low 
masses, where ellipticals are more compact than spirals, disks 
have a large fraction of their mass in gas, so mergers must 
account for dissipation. In a merger of two disks containing 
both gas and stars, the stars are dissipationless and they can- 
not increase their phase space density and so violently relax to 
a distribution with an R e similar to the progenitor disks. Gas, 
on the other hand, can radiate, and tidal torqu es excited during 
a merger can remove its ang ular momentum (Hernquistl l 19891; 
iBarnes & Hernquislll99lLll996l) . The resulting inflows pro- 
duce a dissipational, merger-induced starburst which is com- 
pact (typical size scales ~ 0.5 - 1 kpc). If a significant frac- 
tion of the final stellar mass is formed in this manner, the scale 
length of the stellar component will be much smaller than that 
o f its progenitor . 

lOnorbe et al l (120051) : iRobertson eTall d2006b) and 
iDekel & Coxl ((2006) argued that, because low mass disks are 
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more gas rich (on average) than high mass disks, dissipation 
will be more important in low-mass systems. That is, 
lower-mass ellipticals (the merger products of low-mass 
spirals) should have smaller effective radii relative to their 
progenitor disks, and, since the halo mass distribution is not 
strongly affected by this process, a correspondingly smaller 
dark matter fraction (lower Mdy n /M*) within the stellar 
R e . On the other hand, high-mass disks are observed to be 
gas-poor: at ;§> L* gas fractions become negligible (<C 10%) 
- this is precisely where ellipticals are not more compact than 
disks. Together, iRobertson et al.1 (l2006bl) and iDekel & Coxl 
(2006) argued that the dependence of dissipational fraction 
on mass is sufficient, in principle, to explain the tilt of the FP, 
the size-mass and velocity-dispersion mass correlations of 
ellipticals. 

The importance of gas dynamics and triggered star for- 
mation in mergers is reinforced by observati ons of ul- 
tralum inous infrared galaxies (ULIRGs) (e.g. ISoifer et all 
Il984alfbh . which are invariably associated with mergers in the 
local Universe dJose ph & Wright 1985; Sanders & Mirabel 
1996). The infrared emission from ULIRGs is thought 
to be powered by intense starbursts in their nuclei, origi- 
nating in central concentrations of gas (e.g. IScoville et al.1 
Il986t ISargentet .alj 1 19871 1 19891). whi ch will leave dense 
stella r remnants dKor mendv & Sanders 1992; lHibbard & Yur& 
Il999t iRothberg & Joseph! 120041) , as predicted theoreti- 
cally (Mi hos & Hernquisl Il994ah . Moreover, observations 
of merging systems and gas-rich merger remnants (e.g.. 
Lake & Dressier! I l986t iDoyon et al.1 1 19941: IShier & Fischen 
1998UJames et al.ll999h . as well as post-starburst (E+A/K+A) 
galaxies dGotol 120051) . have shown that their kinematic and 
photometric properties are consistent with them eventually 
evolving into typical ~ L* elliptical galaxies. The correlations 
obeyed by these mergers and r emnants (e.g., iGenzel et al.l 
l200U rRothberg & Joseph 2006a b, and references above) are 
simila r to e.g. the observed fundamental plane and lKormendyl 
(1 19771) relations for relaxed ellipticals, and consistent with 
evolution onto these relations as their stellar populations 
age, as well as the cl ustering and mass density of ellipticals 
dHopkins et al.ll2007bl) . 

Unfortunately, the consequences of the models are less 
clear; as such, there has been some ambiguity regarding 
whether or not recent observations of the FP support or 
disagree with theory. In particular, while merger rem- 
nants may fall on the FP, it is not obvious that a differ- 
ential role of dissipation is in fact responsible for the FP 
tilt in the manner predicted by Robertson et al. (20063), or 
that this applies to all ellip ticals. However, there is hope: 
Mihos & Hernq uisl (Il994al) predicted that these dissipational 
starbursts should leave an observable signature in the surface 
brightness profiles of remnants, in the form of a steep depar- 
ture from the outer Ide Vaucouleursl (1 19481) r'/ 4 -law distribu- 
tion in the inner regions: i.e. a central "extra light" com- 
ponent above the inwards extrapolation of the outer profile. 
Observations have now uncovered distinctive e vidence for 
this two-component structure in local ellipticals (|Kormendvl 
119991; iKormendv etall l2008t jFerrarese et al.1 120061) . classi- 
cal bulges (iBalcells et al.ll2007[). and recent merg er remnants 
dHibbard & Yunlll999t IRothberg & Joseph! 12004 . With the 
combination of HST and ground-based photo metry, it now ap- 
pears that such components are ubiquitous dKormendv et al.1 
2008), with mass ranges and spatial extents comparable to 
those expected from observations of ongoing merger-induced 



starbursts and numerica l simulations. 

In a series of papers, Hopkinsetal] |2008d a.g) (hereafter 
Paper I, Paper II and Paper III, respectively), we used these 
simulations and data to develop and test a method to empiri- 
cally determine the degree of dissipation involved in the for- 
mation of a particular elliptical galaxy - i.e. the mass fraction 
in the stellar remnant of a central, compact nuclear (dissipa- 
tional) starburst. In Paper I, we demonstrated that observed 
merger remnants can be robustly decomposed into two com- 
ponents: an outer, dissipationless (violently relaxed) compo- 
nent with a Sersic law-like profile, comprising the pre-merger 
stars, and an inner, compact starburst remnant, produced in a 
starburst. Combining large ensembles of observations with a 
library of simulations that enabled us to calibrate various em- 
pirical methods, we developed a purely empirical technique 
to separate the inner "excess" owing to the true physical star- 
burst component in the observed surface brightness profile 
from the outer profile. 

In Paper II, we showed that this method - given photom- 
etry of sufficient quality and covering a large dynamic range 
(from < lOOpc to > 20-50kpc) - could be extended to ob- 
served "cusp" ellipticals (i.e. ellipticals with steep nuclear 
profiles), co mmonly believed to be the direct remnants of gas- 
rich mergers dFaber et al.lll997l) . Separating the observed sur- 
face brightness profile in this manner, we demonstrated that 
simulations and independent observations (e.g. distinctions in 
stellar populations evident in kinematics, colors, stellar ages, 
metallicities, or abundances) confirm that the component of 
the elliptical formed via dissipation (in a nuclear starburst) 
could be reliably (statistically) determined. 

In Paper III, we showed that the same methods can robustly 
recover the dissipational starburst remnants in "core" ellip- 
ticals (ellipticals with shallow nuclear profiles). In general, 
even if other processes such as e.g. scattering of stars by a 
binary black hole create the core, their impact on the overall 
starburst component is negligible (by both mass and radius, 
the scales of the core, typically < 30-50pc, are much smaller 
than the mass and size of the starburst). We also showed that 
even if core or other ellipticals have subsequently been mod- 
ified by spheroid-spheroid "dry" re-mergers, profile shape is 
preserved to a sufficient degree that the original nuclear excess 
(i.e. the indicator of the degree of dissipation in the original, 
spheroid-forming merger) remains. 

By applying our methodology to observations of ellipticals 
over a wide range in mass and size, we can, for the first time, 
empirically compare the degree of dissipation (starburst or 
dissipational mass fraction) in ellipticals to their global prop- 
erties and locations on the FP. If a differential effect of dis- 
sipation as a function of mass is the explanation for the FP 
and elliptical scaling relations, as the models suggest, then 
we should be able to see and quantify the signatures of this 
directly in their observed profiles. Therefore, in this paper 
we present a critical examination of the relationship between 
spheroid properties and FP correlations and "extra light" com- 
ponents, in both simulations and observed galaxies. The com- 
bination of a large number of observations, together with an 
ensemble of hydrodynamic gas-rich merger simulations sam- 
pling the entire observed range in e.g. mass, gas content, and 
other properties, enables us to develop and apply new, de- 
tailed, empirical tests of these models for the origin of the FP 
correlations. 

In §|2] we summarize our library of merger simulations, and 
in §[3]we describe the compilation of observations used to test 
the models. In §[4] we review different approaches to fit the 
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surface density profile and recover the physically distinct (dis- 
sipational versus dissipationless) components in merger rem- 
nants. We use a set of simulations to infer how galaxy proper- 
ties are predicted to scale with dissipation in § [5] We then 
compare with observed systems: examining how observed 
sizes and masses scale with gas content in §|6] and the scaling 
relations obeyed at fixed dissipational content in §|7] We com- 
bine the observed dependences on dissipation and gas content 
in § |8] to determine whether this is sufficient to explain the 
tilt and scatter of the FP and its projected correlations. In §|9] 
we compare dynamical and total masses and consider possi- 
ble non-homology effects. In §[10] we examine the impact of 
subsequent re-mergers on the FP correlations. Finally, in § [TT1 
we discuss our results and outline future explorations of these 
correlations. 

Throughout, we adopt a Om = 0.3, ft\ = 0.7, Hq = 
70kms _1 Mpc -1 cosmology, and normalize all observations 
and models accordingly. We note that this h as little affect on 
our conclusions, however. We also adopt a Cha brierl (120031) 
stellar initial mass function (IMF), and convert all stellar 
masses and mass-to-light ratios to this choice. The exact form 
of the IMF systematically shifts the normalization of stellar 
masses herein, but does not substantially influence our com- 
parisons. All magnitudes are in the Vega system, unless oth- 
erwise specified. 

2. THE SIMULATIONS 

The merger simulations we analyze in this paper were 
performed with the parallel TreeSPH code GADGET-2 
dSpringelll2005l). based on a fully conservative formulation 
dSpringel & Her nquist 2002) of smoothed particle hydrody- 
namics (SPH), which conserves energy and entropy simul- 
taneously ewn_when^ooflih2S^SS^_2Y l ve adaptively 
(see e.g.. lHernquistlll993l: lO'Shea et alll2005l) . The simula- 
tions account for radiati ve cooling and optional heating by 
a UV background (as in iKatz et all 119961: iDave et all 119991 
although it is not important for the masses of interest here), 
and incorporate a sub-resolution model of a multiphase in- 
terstellar medium ( ISM) to describe star forma tion and su- 
pernova feedback ( Springel & Hernquist 20031). Feedback 
from supernovae is captured in this sub-resolution model 
through an effective equation of state for star-forming gas, en- 
abling us to stably evolve disk s with arbitrary gas fractions 
(see, e.g. ISprin gel et al. 2005] ISpringel & Hernquistl 120051: 
Rober tson et alll2006allcT) . This is described by the parame- 
ter g'eos, which ranges from q eos = for an isothermal gas with 
effective temperature of 10 4 K, to q eoti = 1 for our full multi- 
phase model with an effective temperature ~ 10 s K. We also 
compare with a subset of simulations whi ch adopt the star for- 
mation an d feed back prescriptions from Mi hos & HernquisH 
dl994allcl 11996b . in which the ISM is treated as a single- 
phase isothermal medium and feedback energy is deposited 
in a purely kinetic radial impulse (for details, see, e.g. 
iMihos & Hernquistl d!994bl) ). 

Although we find that they make little difference to the ex- 
tra light component, most of the simulations include super- 
massive black holes at the centers of both progenitor galaxies. 
The black holes are represented by "sink" particles that ac- 
crete gas at a rate M estimated from the local gas density and 
sound speed using an Eddington-limited prescription based on 
Bondi-Hoyle-Lyttleton accretion theory. The bolometric lu- 
minosity of the black hole is taken to be Lboi = e f Mc 2 , where 
e r = 0.1 is the radiative efficiency. We assume that a small 
fraction (typically w 5%) of Lboi couples dynamically to the 



surrounding gas, and that this feedback is injected into the 
gas as thermal energy, weighted by the SPH smoothing ker- 
ne l This fraction is a fre e parameter, which we determine as 
in iDi Matteo et all d2005l) by matching the observed M B h - o 
relation. For now, we do not resolve the small-scale dynam- 
ics of the gas in the immediate vicinity of the black hole, but 
assume that the time-averaged accretion rate can be estimated 
from the gas properties on the scale of our spatial resolution 
(roughly ks 20 pc, in the best cases). In any case, repeating 
our analysis for simulations with no black holes yields identi- 
cal conclusions. 

The proge n itor galaxy models are described in 
Springel et al. (2005), and we review their properties 
here. For each simulation, we generate two stable, isolated 
disk galax i es, ea ch with an extended dark matter halo with a 
iHernquisll d 19901) profile, motivated by cosm ological simula- 
tions dNavarro et al]|1996l : [Busha et al.l l2005). an exponential 
disk of gas and stars, and (optionally) a bulge. The galaxies 
have total masses M vir = V^J{\0GH[z\) for an initial redshift 
Z, with the baryonic disk having a mass fraction m& = 0.041, 
the bulge (when present) having m\, = 0.0136, and the rest of 
the mass in dark matter. The dark m atter halos are as s igned a 
concentration parameter scaled as in lRobertson et alj (2006c) 
appropriately for th e galaxy mass and redshift following 
Bullo cket all d2001l) . We have also varied the concentration 
in a subset of simulations, and find it has little effect on 
our conclusions because the central regions of the galaxy 
are baryon-dominated. The disk scale-length is computed 
based on an assumed spin parameter A = 0.033, chosen to be 
near the mode in the A distribution measured in simulations 
dVitvitska et alll2002l) . and the scale-length of the bulge is set 
to 0.2 times this. Modulo explicit variation in these parame- 
ters, these choices ensure that the initial disks are consistent 
with e.g. the observed baryonic Tully-Fish er relation and 
estimated halo-galaxy mass scaling laws (Bell & d e Jond 
1200 U iKormendv & Freemanll2004t iMandelbaum et alj|2006l 
and references therein). 

Typically, each galaxy initially consists of 168000 dark 
matter halo particles, 8000 bulge particles (when present), 
40000 gas and 40000 stellar disk particles, and one black hole 
(BH) particle. We vary the numerical resolution, with many 
simulations using twice, and a subset up to 128 times, as many 
particles. We choose the initial seed mass of the black hole 
either in accord with the observed Mbh-c relation or to be 
sufficiently small that its presence will not have an immediate 
dynamical effect, but we have varied the seed mass to iden- 
tify any systematic dependencies. Given the particle num- 
bers employed, the dark matter, gas, and star particles are all 
of roughly equal mass, and central cusps in the dark matter 
and bulge are reasonably well resolved. The typical gravita- 
tional softening in the simulations is ~ 20-50pc in the < L* 
systems of particular interest here, with a somewhat higher 
~ 50- lOOpc in the most massive systems (yielding an effec- 
tively constant resolution ~ 0.01 R e in terms of the effective 
radius). In Paper I and Paper II we demonstrate that this is suf- 
ficient to properly resolve not only the mass fractions but also 
the spatial extent of the extra light components of interest here 
(resolution may become an issue when attempting to model 
the very smallest galaxies, with R e < 100 pc and L < 0.01 L*, 
but this is well below the range of the observations of interest 
here). The hydrodynamic gas smoothing length in the peak 
starburst phases of interest is always smaller than this gravi- 
tational softening. 

We consider a series of several hundred simulations of col- 
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liding galaxies, d escribed in iRobertson et al.l d2006bllcl) and 
ICox et al.l(l2006 a b). We vary the numerical resolution, the or- 
bit of the encounter (disk inclinations, pericenter separation), 
the masses and structural properties of the merging galax- 
ies, initial gas fractions, halo concentrations, the parameters 
describing star formation and feedback from supernovae and 
black hole growth, and initial black hole masses. 

The progenitor galaxies have virial velocities V vu = 
55,80,113,160,226,320, and 500kms-', and redshifts z = 
0,2, 3, and 6, and the simulations span a range in final 
spheroid stellar mass M» ~ 10 s - 10 I3 M Q , covering essen- 
tially the entire range of the observations we consider at all 
redshifts, and allowing us to identify any systematic depen- 
dencies in the models. We consider initial disk gas fractions 
by mass of / gas = 0.05, 0.1, 0.2, 0.4, 0.6, 0.8, and 1.0 (de- 
fined as the fraction of disk baryonic mass which is gas) for 
several choices of virial velocities, redshifts, and ISM equa- 
tions of state. The results described in this paper are based 
primarily on simulations of equal-mass mergers; however, by 
examining a small set of simulations of unequal mass merg- 
ers, we find that the behavior does not change dramatically for 
mass ratios to about 3:1 or 4:1. The mass ratios we study are 
appropriate for the observations of ellipticals used in this pa- 
per, which are only formed in our simulations in major merger 
events. At higher mass ratios, the resu lt is a small bulge in 
a still disk-dominate d galaxy (see e.g. Younger et all 120081 
Hopkins et al. 2008c f), which we do not study here (although 
in general our conclusions should still apply, so long as the 
bulges of interest are "classical" bulges formed in mergers). 

We also briefly consider in § [10] a subset of spheroid- 
spheroid "re-mergers," representative of gas-poor or "dry" 
spheroid-spheroid mergers of elliptical galaxies. In these 
mergers, we collide two remnants of previous disk-disk merg- 
ers, in order to explore how their properties are modified 
through re-merging. We typically merge two identical rem- 
nants (i.e. two identical copies of the remnant of a given disk- 
disk merger), but have also explored re-mergers of various 
mass ratios (from 1:1 to « 4:1), and mixed morphology re- 
mergers (i.e. merging an elliptical remnant with an un-merged 
gas-rich disk). In the former case, we generally find a similar 
division in mass ratio at which a major merger is significant. 
In the latter, we find that the properties are more akin to those 
of other gas-rich (disk-disk) mergers, and the remnant should 
for most purposes should still be considered the direct prod- 
uct of a gas-rich merger. In the re-merger series, we vary the 
orbital parameters, both of the initial gas-rich merger and re- 
merger, and consider systems with a similar range of initial 
gas fractions in the progenitor disks (of the original gas-rich 
merger), / gas = 0.05, 0.2, 0.4, 0.8. The re-mergers span a 
similar range in virial velocities and final stellar masses to the 
gas-rich mergers. 

Each simulation is evolved until the merger is complete and 
the remnants are fully relaxed, typically ~ 1 -2 Gyr after the 
final merger and coalesc ence of the BHs. We then analyze 
the remnants following ICox et all (l2006bl) . in a manner de- 
signed to mirror the methods typically used by observers. For 
each remnant, we project the stars onto a plane as if observed 
from a particular direction, and consider 100 viewing angles 
to each remnant, which uniformly sample the unit sphere. 
Given the projected stellar mass distribution, we calculate the 
iso-density contours and fit ellipses to each (fitting major and 
minor axis radii and hence ellipticity at each iso-density con- 
tour), moving concentrically from r = until the entire stellar 
mass has been enclosed. This is designed to mimic obser- 



vational iso photal fitting algorithms (e.g. lBender et al.ll 1 987b 
lBenderlfl988h . The radial deviations of the iso-density con- 
tours from the fitted ellipses are expanded in a Fourier series 
in the standard fashion to determine the boxyness or diskyness 
of each contour (the ct\ parameter). Throughout, we show pro- 
files and quote our results in terms of the major axis radius. 
For further details, we refer to lCox et"aT1 d2006bl) . 

We directly extract the effective radius R e as the projected 
half-mass stellar effective radius, and the velocity dispersion 
a as the average one-dimensional velocity dispersion within 
a circular aperture of radius R e . This differs from what is 
sometimes adopted in the literature, where R e is determined 
from the best-fitting Sersic profile, but because we are fitting 
Sersic profiles to the observed systems we usually quote both 
the true effective radius of the galaxy and effective radii of 
the fitted Sersic components. Throughout, the stellar mass M* 
refers to the total stellar mass of the galaxy, and the dynamical 
mass Mdy n refers to the traditional dynamical mass estimator 

M dyn = k-^, (4) 

where we adopt k = 3.8 (roughly what is expected for a 
iHernquisll (1990) profile, and the choice that most accurately 
matches the true enclosed stellar plus dark matter mass within 
R e in the simulations; although this choice is irrelevant as 
long as we apply it uniformly to both observations and simu- 
lations). When we plot quantities such as R e , a, and Md yn , 
we typically show just the median value for each simula- 
tion across all ~ 100 sightlines. The sightline-to-sightline 
variation in these quantities is typically smaller than the 
simulation-to-simulation scatter, but we explicitly note where 
it is large. 

3. THE DATA 

We compare the simulations to and test our predictions 
on an ensemble of observed surface brightness profiles of 
ellipticals, described in Paper II and Paper III. Specifi- 
cally, we consider three samples of ellipticals and a com- 
pilation of remnants of recent gas-rich me rgers. The first 
is the V-band Virgo elliptical survey of Kor mendv et alj 
(2008), based on the complete sampl e of Virgo galaxies 
down to extremely faint s ystems in iBinggeli et al.l (1985) 
(the s a me sample studied in | C6te et al.ll2006l; iFerrarese et alj 
2006). Korme ndvetail (120081) ' combine observations from a 
large number of sources (including Bender et al. 1988, 2007; 
I Caon et al.l 1 1 9901 |1994t iDavis et alJ 119851; iKormendy et alj 
20051: lLaueri | 1985[ lLauer et al.lll995Ll2005l iLiu et al.ll2005F 
iPeletier et al.ll990» and new photometry from McDonald Ob- 
servatory, the HST archive, and the SDSS for each of their ob- 
jects which (after careful conversion to a single photometric 
standard) enables accurate surface brightness measurements 
over a wide dynamic range (with an estimated zero-point ac- 
curacy of ±0.04 V mag arcsec" 2 ). Typically, the galaxies in 
this sample have profiles spanning ~ 12-15 magnitudes in 
surface brightness, corresponding to a range of nearly four or- 
ders of magnitude in physical radii from ~ lOpc to ~ 100 kpc, 
permitting the best simultaneous constraints on the shapes of 
both the outer and inner profiles of any of the objects we study. 
The profiles include e.g. ellipticity, a\ja, and g— z colors as 
a function of radius. Unfortunately, since this is restricted to 
Virgo ellipticals, the number of galaxies is limited, especially 
at the intermediate and high end of the mass function. 

We ther e fore a dd surface brightness p rofiles from 
lLauer et all d2007al) . further supplemented by iBender et alj 
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(1988). lLauer et all d2007al) compile V-band measurements 
of a large number of nearby systems for which HST imag- 
ing of the galacti c nuclei is available. These include the 
lLauer et alj d2005l) WFPC2 data-set, the lLaine et all (120031) 
WFPC2 BCG sample (in which the objec ts are specifically 
selected as brigh test cluster g alaxies from Postman & Lauer 
(1995)), and the lLauer et alld 19951) and iFaber etalj d 19971) 
WFPC1 compilations. Details of the treatment of the profiles 
and co nversion to a single standard are given in Lau er et al.1 
d2007al) . The sample includes ellipticals over a wide range of 
luminosities, down to Mb ~ -15, but is dominated by interme- 
diate and giant ellipticals, with typical magnitudes Mb < — 18. 
This therefore greatly extends our sampling of the intermedi- 
ate and high-mass end of the mass function, but at the cost 
of some dynamic range in the data. The HST images alone, 
while providing information on the central regions, typically 
extend to only ~ 1 kpc outer radii, wh ich is insufficient to 
fit the outer profile. lLauer et ail (l2007al) therefore combine 
these data with ground-based measu rements from a number 
of sources (see the references for the iKormendv e t al. (2008) 
sample) to construct profiles that typically span physical radii 
from ~ lOpc to ~ 10 — 20 kpc. Althou gh the composite pro- 
files were used in lLauer et all d2007al) to estimate effective 
radii, they were not actually shown in that paper. It should 
also be noted that there is no single criterion that characterizes 
galaxies included in this sample, but they generally comprise 
luminous nearby ellipticals and SO galaxies for which detailed 
imaging is available. We emphasize that issues of complete- 
ness and e.g. environment are not important for any of our 
conclusions. 

We o ccasionally supplement the profi les from lLauer et al.1 
d2007al) wit h additional profiles used in iBenderet all d 19881 
I1992H19931 [T994). and in some cases subsequently updated. 
These are more limited: typically the profiles cover ~ 7 mag- 
nitudes in surface brightness, extending from ~ 30-50pc out 
to ~ a few kpc (typically ~ 3 kpc in low-luminosity systems, 
and ~ 15 kpc in the brightest systems, sufficient for accept- 
able, but not strong constraints on the outer profile shapes). 
However, the measurements are usually in each of the V, R, 
and / bands, and hence allow us to construct multicolor sur- 
face brightness, ellipticity, and 04/0 profiles. We use this to 
estimate e.g. the sensitivity of the fitted parameters and galaxy 
profiles on the observed waveband and on the quality and dy- 
namic range of the photometry. 

We also consider t he sample of local remnan ts of recent 
gas-rich mergers from Rothberg & Joseph (2004) with which 
we compare our simulati ons in Paper I. For these objects, 
Rothberg & Joseph ( 2004) compile K-bund imaging, surface 
brightness, ellipticity, and a^ja profiles, where the profiles 
typically range from ~ lOOpc to ~ 10-20 kpc. These span a 
moderate range in luminosity (including objects from Mk ~ 
-20 to M K ~ -27, but with most from M K ~ -24 to M K ~ -26) 
and a wide range in merger stage, from ULIRGs and (a few) 
unrelaxed systems t o she ll ellipticals. As demonstrated in 
Rothberg & Joseph! d2004l) and argued in Paper I, these sys- 
tems will almost all become (or already are, depending on 
the classification scheme used) typical ~ L» ellipticals, with 
appropriate phase space densities, surface brightness profiles, 
fundamental plane relations, kinematics, and other properties. 
For a detailed discussion of the modeling of these systems 
and the profiles themselves, we refer to Paper I (all of the re- 
sults shown for these systems are derived therein). We show 
the results from Paper I here in order to test whether observed 
merger remnants (not just our simulations) obey the same cor- 



relations we study in ellipticals. In particular, this allows us to 
provide an additional empirical check of the simulations and 
the argued continuity of merger and elliptical populations. 

Because the generally accepted belief is that core ellipti- 
cals were not directly formed in gas-rich major mergers but 
were subsequently modified by dry re-mergers, we have re- 
peated our analysis considering just those ellipticals which are 
confirmed via HST observations as being either cusp or core 
ellipticals, and designate the two populations as two separate 
sub-samples throughout. However, as we discuss below, treat- 
ing them as a single population (or including all ellipticals, 
even those without HST observations, in our sample) makes 
no difference to our conclusions. We include all the confirmed 
gas-rich merger remnants, but note there are a small number 
of extreme unrelaxed cases for which sharp features in the 
surface brightness profiles prevented derivation of meaning- 
ful quantities (note, however, as shown in Paper I, that almost 
all of the objects in this sample are sufficiently well-relaxed at 
the radii of interest for our fitting). We exclude spheroidals, as 
they are not believed to form in major mergers as are ellipti- 
cals (e.g. lKormendvlll985L[T987t IKormendv et al.ll2008h , and 
in any case they dominate at extremely low masses where our 
simulations do not sample the population (they also predomi- 
nate as satellite galaxies, whose effects we do not model). 

We also exclude SO galaxies. This is not because of a phys- 
ical distinction: observations suggest that these likely form 
a continuous family with the low-luminosity cusp ellipticals, 
and in fact a number of our simulated gas-rich merger rem- 
nants would, from certain viewing angles, be classified as 
SOs (although we exclude disk-dominated simulation rem- 
nants from our comparisons here). However, in order to derive 
e.g. the parameters of the outer, violently relaxed profile and 
central extra light, it would be necessary to remove the con- 
tribution of the large-scale disk from the surface brightness 
profiles of these objects. Our two-component (outer dissipa- 
tionless and inner dissipational) Sersic models (described in 
§ @|i then become three-component fits, and the degeneracies 
involved with three independent components, even with our 
best data and simulations, are so large as to render the results 
meaningless. We have, however, re-visited all of the SOs in 
these samples in light of our results, and find that they are, in 
all cases, consistent with our predicted and observed trends. 
However, it is too difficult to infer these trends directly from 
the SOs themselves without ideal disk subtraction. 

This yields a final sample of rs 180 unique elliptical galax- 
ies, and k, 50 confirmed remnants of gas-rich mergers. Most 
of the sample spans a range of three orders of magnitude in 
stellar mass, from < 0.1M, to ~ 10M*, and a wide range in 
extra light properties. There is, of course, some overlap in the 
samples that define our compilation; we have ~ 600 surface 
brightness profiles for our collection of unique ellipticals, in- 
cluding (for many objects) repeated measurements in multiple 
bands and with various instruments. This turns out to be quite 
useful, as it provides a means to quantify error estimates in 
fits to these profiles. The variations between fit parameters 
derived from observations in different bands or made using 
different instruments are usually much larger than the formal 
statistical errors in the fits to a single profile. There are no ob- 
vious systematic effects (i.e. systematic changes in profile fits 
from V to / band), but as demonstrated in Paper I the effects 
of using different bands or changing dynamic range (from dif- 
ferent instruments) can be complex, depending on the struc- 
ture and degree of relaxation of the outer regions of a system. 
On the other hand, there are well-relaxed objects for which 
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almost no significant change in the fits occurs from band to 
band. It is therefore useful to have multiple observations of 
the same system, as it allows us to get some idea of how sen- 
sitive our fits are to differences in e.g. the choice of observed 
wavelength or dynamic range from instrument to instrument. 

In Paper II and Paper III, we present the results of our fits to 
each elliptical in our sample; we use these values throughout 
this paper. For sources with multiple independent observa- 
tions, we define error bars for each fit parameter representing 
the ~ 1 a range in parameters derived from various observa- 
tions, typically from three different surface brightness profiles 
but in some cases from as many as w 5 - 6 sources (where 
there are just 2 sources, the "error" is simply the range be- 
tween the two fits). In many cases the different observations 
are comparable; in some there are clearly measurements with 
larger dynamic range and better resolution: the errors derived 
in this manner should in such cases be thought of as the typ- 
ical uncertainties introduced by lower dynamic range or less 
accurate photometry. 

In terms of direct comparison with our simulations, the data 
often cover a dynamic range and have resolution compara- 
ble to our simulations, provided we do not heavily weight the 
very central (< 30pc) regions of HST nuclear profiles. Ex- 
perimenting with different smoothings and imposed dynamic 
range limits, we find it is unlikely that resolution or seeing 
differences will substantially bias our comparisons. They can 
introduce larger scatter, however: the robustness of our results 
increases considerably as the dynamic range of the observed 
profiles is increased. 

We have converted all the observations to physical units 
given our adopted cosmology, and compile global parame- 
ters (where not available in the original papers) including 
e.g. kinematic properties, luminosities, and black hole masses 
from the literature. We determine stellar masses ourselves in a 
uniform manner for all the objects, based on their total A'-band 
lumin osities and (B—V ) color-dependent mass-to-light ratios 
from iBell et al.l (120031) . corrected for our adopted IMF. We 
have also repeated our analysis using stellar masses derived 
from a mean M/L as a function of luminosity or from fitting 
the integrated UBVRIJHK photometry of each object to a sin- 
gle ste llar population with the models of iBruzual & C hariot 
(2003), and find this makes no difference to our conclusions. 

Throughout, we will usually refer interchangeably to the 
observed surface brightness profiles in the given bands and the 
surface stellar mass density profile. Of course, stellar light is 
not exactly the same as stellar mass, but in Paper I and Paper 
II, we consider the differences between the stellar light and the 
stellar mass density profiles as a function of time, wavelength, 
and properties of the merger remnant, and show that the opti- 
cal (V) and Zf-band results introduce little bias (i.e. are good 
tracers of the stellar mass). It is important to note that while 
we are not concerned about the absolute normalization of the 
profile (i.e. mean M/L), since we derive total stellar masses 
separately from the integrated photometry, we must account 
for systematics that might be induced by change in M/L as 
a function of radius. The profiles in optical bands such as V 
require more care when the system is very young (< 1 - 2 Gyr 
after the major merger-induced peak of star formation), and 
there can be considerable bias or uncertainty owing to stellar 
population gradients and dust. However, once the system is 
relaxed, the optical bands become good proxies for the stellar 
mass distribution, with < 20% variation in our M/L over the 
entire fitted range of radii in our simulations, in good agree- 
ment with the simple expectation based on the observed weak 



color gradients in most of the observed systems. 

In fact, in Paper I and Paper II we demonstrate that once the 
system reaches intermediate age, the bias in e.g. B or V band 
is often less than that in K band, because systems tend to be 
both younger and more metal rich in their centers. In /f-band, 
these both increase L/M, leading to a (small) systematic bias. 
Our simulation results indicate that our merger remnant sam- 
ple (observed in /f-band) are, on average, robust in this sense, 
but they should be treated with care, especially in the most ex- 
treme cases (namely the few LIRGs and ULIRGs in the sam- 
ple). In optical bands, however, age and metallicity gradients 
have opposite effects (younger age increases L/M, but higher 
metallicity decreases L/M), and they tend to mostly cancel. 
Since essentially all of our ellipticals observed in the optical 
bands are older than this stellar population age (even in their 
centers), and they have been carefully vetted and either cor- 
rected or (where correction was too di fficult) excluded for the 
effects of e.g. dust lanes in the sources dKormendv et alj2008t 
lLauer et aLll2007al) . we are not concerned that significant bias 
might persist. Furthermore, comparison of systems observed 
in different bands demonstrates that our conclusions are un- 
changed (modulo small systematic offsets) regardless of the 
observed bands in which we analyze these systems. 

4. INFERRING THE DEGREE OF DISSIPATION FROM OBSERVED 
PROFILES 

In ord er to test the theoretical models from Rob ertson et al.l 
(2006b) and others using observed systems, we require some 
empirical means to estimate the amount of dissipation in- 
volved in forming elliptical galaxies. The simulations imply 
that the light profile should be considered to be the sum of two 
physically distinct components. The outer, violently relaxed, 
dissipationless component is made from stars that formed in 
the disks prior to the final stages of a merger, and therefore 
evolve in a collisionless manner. The inner, "dissipational" 
component is comprised of stars produced in a central star- 
burst, from gas which loses its angular momentum in the 
merger. This gas is channeled into the center of the remnant 
and is converted into stars on a short timescale (effectively in 
an almost fixed background potential set by the already nearly 
relaxed outer, dissipationless component). The "degree of dis- 
sipation" is, therefore, effectively measured by the mass frac- 
tion of the central, dissipational starburst component, which 
causes R e to contract and gives rise to the effects described 
above. In Paper I, Paper II, and Paper III, we develop and test 
a methodology for decomposing observed surface brightness 
profiles into these two components. We refer to those papers 
for details and verification that this approach is reliable when 
applied to either simulations or observations, but we briefly 
review the methodology here. We consider two different ob- 
servable proxies for the amount of dissipation. 

First, we directly fit the surface brightness profile to deter- 
mine the mass fraction in a central "extra light" or starburst 
component, which we denote as /extra- We demonstrate in Pa- 
per I that this quantity can be measured by fitting the total 
observed surface brightness profile of merger remnants to the 
following two component model: 

/tot=4xtraexp{-&i(^^)}+/ exp{-&„(^) / j, (5) 

with the inner part (7 extra , r extra ) serving as a proxy for the 
central dissipational component and the outer component (7 , 
r a ) representing the dissipationless outer profile. The "extra 
light," in other words, is the mass fraction in a high surface 
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FIG . 1 . — Left: Surface mass density of a typical merger remnant from our 
simulation library (black), decomposed into stars formed prior to the final 
merger (which are then violently relaxed; red) and stars formed in the dissi- 
pational starburst (blue). The Sersic index fitted to the pre-starburst compo- 
nent alone is shown, with the mass fraction of the starburst component. We 
show two examples, one very gas-rich (top) and one less gas-rich (bottom). 
Right: Two-component (Sersic plus cusp or extra light) fit (inner exponen- 
tial and outer Sersic) to the total light profile, with the Sersic index of the 
outer component and mass fraction of the inner component, and rms scatter 
(Afj,) about the fit. Our two-component, cusp plus Sersic function fit (right) 
accurately recovers the profile of the violently relaxed component and mass 
fraction of the starburst component. 

brightness central component that rises above the extrapola- 
tion of the outer, intermediate n s Sersic law-like profile (which 
in any reasonable model will be dominated by the violently 
relaxed stars). In Paper I and Paper II, we show that this 
form provides, on average, a reliable physical decomposition 
of the mass distribution in remnants of simulated mergers be- 
tween gas-rich galaxies, including mass fractions, effective 
radii, and profile shapes of the two components. Analogous 
decompositions have been applied to e.g. observed bulges in 
disk galaxies dBalcells et al J 12007) and to Virgo ellipticals, 
motivated by purely empirical considerations including sur- 
face brightness profile shapes, isophotal sha pe and kinematic 
profiles, and central s tellar populations (see Kormendvl l 1 9991 : 
iKormendV et al. 2008| and references therein), a s well as re- 
cent merger remnants dRothberg & Josephll2004l) : we refer to 
Paper II for a detailed discussion and comparison with similar 
methodologies from the literature. 

Figure Q] illustrates this procedure for two typical merger 
remnants from in our simulation library. We show the to- 
tal surface brightness profile, as well as the known physical 
decomposition of that profile into a pre-starburst, violently 
relaxed component and a dissipational starburst component 
(here, we use a cut in time about the peak starburst to identify 
the "starburst" stars, but as shown in Paper I, our results are 
not sensitive to the exact choice). We compare this with the 
results obtained by fitting the total profile (temporarily ignor- 
ing the known physical decomposition) to a two-component 
model of the form given by Equation ©. In both cases, the in- 
ner component of our fit is a good proxy for the physical star- 
burst light distribution, accurately capturing its mass fraction, 
effective radius, and shape where it is important; i.e. where 
it contributes significantly to the light profile. Likewise, the 
shape (Sersic index), radius, and mass fraction of the dissi- 
pationless component are accurately recovered by this purely 
empirical method. 

We apply this to observed merger remnants in Paper I and 
cusp ellipticals in Paper II, and show that it yields good fits 
with reliable results. In particular, we find that where there 



is independent evidence for a distinct stellar population (i.e. 
a starburst superimposed on a more extended, older stellar 
population), this method recovers the empirically estimated 
starburst population. We demonstrate that other indicators, 
including e.g. stellar population gradients, kinematics, and 
isophotal shapes all support the physical nature of these de- 
compositions. 

In Paper III, we apply this methodology to remnants of gas- 
poor, spheroid-spheroid re-mergers, and show that, while the 
profiles are somewhat smoothed out, the general profile shape 
is preserved, and applying this procedure will still reliably re- 
cover the original breakdown between dissipational and dis- 
sipationless components (i.e. the starburst component formed 
in the original, spheroid-forming merger). We extend our 
analysis to a large sample of "core" ellipticals, and similarly 
demonstrate the accuracy of this decomposition and agree- 
ment with stellar population models (where available); even 
if "scouring" by a binary black hole has flattened the profiles 
on < 30-50pc scales, the bulk of the "extra light" has not 
been strongly affected and can still be recovered. 

This emphasizes an important caveat: the merger history 
and series of induced dissipational events in a given galaxy 
may be more complex t han a single or couple of idealize d 
major mergers (see e.g. Kobavashi 2004; Naab et alj|2.007h . 
Moreover, merger-induced starbursts may not be the only 
source of dissipation; for example, stellar mass loss may re- 
plenish the gas supply and le ad to new dissipational bursts 
(see e.g. lCiotti & Ost riker 2007). For our purposes, however, 
all dissipational star formation will appear similar when ob- 
served and have the same effects - it is convenient to simulate 
idealized cases, but our results should most appropriately be 
considered a measurement of the integrated amount of dissi- 
pation (regardless of other details of the merger and dissipa- 
tional history). Experiments with e.g. more complex merger 
histories, simulations of multiple simultaneous mergers, and 
series of dissipational events after "rejuvenation" suggest that 
in this sense, our results are robust and independent of the de- 
tailed history. Furthermore, the agreement between these esti- 
mators and independent observational constraints from stellar 
populations, kinematics, and isophotal shapes suggest that the 
recovery is robust. 

Second, we consider the physical starburst mass fraction / s b 
implied by directly fitting simulated surface brightness pro- 
files to the observations. In detail, we demonstrate in Paper I, 
Paper II, and Paper III that in our large library of simulations 
(with varying masses, gas fractions, orbital parameters, stellar 
and black hole feedback prescriptions, and other properties) 
there are remnant surface brightness profiles which (modulo 
small offsets in the exact normalization) agree well with the 
observed profiles over the entire dynamic range of the data. 
Considering the best fit simulation in each case, we find that 
there are almost invariably simulations with similar profiles 
- variance less than Afj, = O.lmagarcsec" 2 of the simulated 
profile with respect to the observed profile (comparable to the 
inherent point-to-point scatter obtained with arbitrary spline 
fits to the simulated or observed profiles). 

These matches are insensitive to e.g. orbital parameters, 
disk initial conditions, or prescriptions for feedback and star 
formation (as expected given the relative independence of the 
inferred dissipational fractions on details of a given merger 
history), but good fits to a given observed profile are obtained 
for only a narrow range in the physical starburst mass frac- 
tion / s b (what dependence there is on other parameters tends 
to indirectly reflect this - varying orbital parameters, for ex- 
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FIG. 2. — Surface density profiles of two simulated gas-rich merger rem- 
nants (the same from Figure [T}, as a function of radius r (bottom) or r'/ 4 
(top). We show both the stellar' profile and the dark matter (halo) profile, and 
label the effective radii of each (vertical lines). Solid line shows a case with 
a moderate degree of dissipation - a mass fraction in the final, centrally con- 
centrated and dissipational merger-induced starburst / s t, = 0.08. Dotted line 
shows a very gas-rich merger remnant, with / s b = 0.31. The profile shapes are 
similar (the non-homology effects are weak), and the halos are nearly identi- 
cal, but the effective radii of the stellar distributions are quite different, owing 
to the dense central concentration of mass from the starburst in the latter case. 



ample, can alter the time required for a merger and therefore 
the mass of gas still available at the time of the final merger to 
participate in the central starburst). We demonstrate that we 
can define a robust and tightly constrained best-fit starburst 
mass fraction, in a \ 2 sense, from fitting each observed pro- 
file to the entire set of simulations. These estimates agree well 
with the "extra light" mass fractions / extr a directly determined 
by profile fitting, and with the other indicators we consider, 
lending further confidence to the view that / extra does indeed 
represent a physically meaningful indicator of dissipation. 

5. THE IMPACT OF DISSIPATION: A CASE STUDY 

In previous pa pers (jRobertson et al.l l2006bt iDekel & Coxl 
2006) argued and Hopki ns et al.l d2008dlah demonstrated that 
increasing the quantity of gas available during the final stages 
of a merger (i.e. increasing the dissipative component of the 
merger which can collapse to small scales and form stars in 
a central starburst) leads to more compact remnants. Because 
of its importance for our analysis here, we use a subset of 
our simulations to highlight the physical significance of this 
result. 

Figure [2] shows two illustrative surface density profiles of 
major merger remnants - the same cases shown in Figure [T] 
The initial galaxies were otherwise identical ~ L* (Milky 
Way-like) spirals placed on a random orbit, except that in 



the first example, the progenitor disks had initial gas fractions 
(~ 1.5 Gyr before the final merger) of 20%, and in the other 
case, of 80%. In the two panels, we show the stellar mass 
density profile of the final (relaxed) merger remnant, and sep- 
ara tely decompose this i nto two physical components follow- 
ing Mihos & H ernquistl (Il994ah : namely an outer component 
comprised of those stars that formed prior to the final coa- 
lescence of the two disks, and were therefore violently re- 
laxed in the merger, and those stars which formed in the (by 
then roughly static) center of the galaxy in the final, merger- 
induced starburst. They correspond to ~ 8% and ~ 30% of the 
final total stellar mass in the two simulations, respectively. 

In both cases, the violently relaxed stars produce an ex- 
tended profile that is well-described by a Sersic law with in- 
dex n s ~ 4. Unsurprisingly, in the example with an initial 
gas fraction ~ 80%, the central light excess reflecting the fi- 
nal starburst component is much more massive - however, the 
profile shape is similar, and this has relatively little effect on 
the outer light profile (it is slightly more compact, because the 
dense central mass concentration means stars scattered to an 
orbit of a given energy will remain at slightly smaller radii in 
this potential, but this effect is small enough that it is not im- 
portant). It is clear from Figure Q] that the differences in the 
two profiles primarily reflect the mass fraction of the central, 
dissipational component, rather than any change in the dissi- 
pationless component. We also show the surface mass density 
of the final dark matter halo. As expected, baryons dominate 
the density within the central ~ 1 — 2 kpc of the system, but 
then trail off. Clearly, the halo mass profile is not very differ- 
ent in the two cases. 

Despite the fact that the central mass concentration does 
not strongly influence the shape and size of the dissipation- 
less remnants of the merger (both the pre-merger, violently 
relaxed stars which constitute the dissipationless stellar com- 
ponent of the merger, and the dark matter halo), it nevertheless 
represents a non-negligible component of the stellar mass, and 
therefore the effective radius is smaller in the case with 80% 
gas (~ 30% of its mass in the final, merger-induced compact 
starburst component). The resulting difference in R e (for such 
a large difference in gas fraction) is dramatic: the effective ra- 
dius in the case with / s b = 0.08 is « 6.8 kpc, and for the case 
with / s b = 0.31 it is » 1.4 kpc. In the latter example, nearly 
half the mass is actually in the dissipational component, so 
the effective radius becomes small. 

How does this effect relate to the ratio of dark matter to 
stellar mass? Figure [3] plots the cumulative dark matter mass 
fraction as a function of radius, 

= M DM (< r) 
M DM «r)+M*«r) 
(for convenience we ignore the negligible mass of gas which 
survives the merger). We also plot the (trivially related) ratio 
of total enclosed mass to total enclosed stellar mass, 
M t Qt = M DM (<r)+M*(<r) 
M* ~ M*(<r) 
It is clear that the difference as a function of radius, at least 
from ~ lOOpc to > 50 kpc, is not large between the two simu- 
lations (which should not be surprising: the difference in their 
stellar mass density profiles is primarily evident at small r, 
where a relatively small fraction of the total mass is enclosed). 
However, if we take the value of these ratios within the effec- 
tive radii R e of the stellar light, we obtain two quite different 
results. In the gas-rich case, the large central light concentra- 
tion yields a small effective radius for the stellar light, within 
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FIG. 3. — Top: Ratio of dark matter mass to total mass enclosed within the 
radius of the elliptical isophotal aperture with major-axis radius r, for the two 
simulations shown in Figureff] Vertical lines denote the effective radii of each 
simulation. Middle: The corresponding ratio of dark matter mass to stellar 
mass enclosed. Bottom: The ratio of the radius-dependent dynamical mass 
estimator at r, 3.8r<r(r)-/G (where o(r) is measured in a narrow annulus at 
r) to the stellar mass enclosed. There is almost no difference in the scaling of 
the two simulations - i.e. the systems appear (in this sense) homologous, but 
the change in R e is such that evaluating these quantities at R e will give rather 
different answers for the two systems. 

which, for both remnants, the baryons dominate the mass. In 
the case with less dissipation, and therefore a weaker cen- 
tral component and larger effective radius, the ratio is mea- 
sured further out (at the larger R e ), where in both simulations 
there are about equal masses enclosed of both dark matter and 
baryons. 

To make a better analogy with observations, we can con- 
struct a dynamical mass estimator 



M dyn «r) = ka(r) 2 r/G, 



(8) 



where we adopt the factor of k = 3.8 for convenience, as noted 
in § E] (making the estimator not so far from the true mass), 
and a(r) is the mean (light-weighted) projected velocity dis- 
persion within a narrow annulus at r. We could also use the 
mean (light-weighted) dispersion within r, (a(< r)), which 
gives a similar result (but traces the enclosed mass less faith- 
fully). Figure [3] compares the ratio Ma yn /M* as a function of 
radius, which yields a similar (albeit noisier) result as con- 
sidering Mtot/M*. Again, it is primarily the difference in the 
effective radii at which this quantity is evaluated, rather than 
the absolute value at fixed r, that drives the difference between 
the gas rich and gas poor mergers. 

We obtain the same result using either the true total en- 
closed mass or a dynamical mass estimator. This implies that 
there is not substantial kinematic or "traditional" structural 
non-homology between our two merger remnants. We can, 




r [kpc] 



FIG. 4. — Ratio of the radius-dependent dynamical mass estimator at r, 
3.&rcr(r) 2 /G (where <r(r) is measured in a narrow annulus at r) to the true 
total (stellar plus dark matter) mass enclosed (M tm e M), for the simulations 
in Figure f5] (shown in the same style; solid and dotted denote low and high 
degrees of dissipation). Vertical lines show the effective radii of each simula- 
tion. The amount of dissipation makes almost no difference in M^ yn /Mi ms: at 
any radius, nor is there a significant difference when evaluated at the effective 
radii of either system. The "homology assumption" that Mi yn oc Mt mc is valid 
in these simulations. 

however, check this directly. 

We compare the ratio M^ yn (r)/M lot (r) as a function of 
radius in Figure [4] The quantity Md yn provides a good 
tracer of M tot , their ratio changing by less than a factor of 
~ 2 over three orders of magnitude in radius, and the ra- 
tio varies only slightly between the two simulations. The 
overall structure of the system has not been significantly 
modified - and in particular over the range of the effective 
radii of the two mergers, there is almost no dependence of 
Md yn (r)/M tot (r) on either the simulation gas content or ra- 
dius, consistent with "isothermality" implied by observational 
constraints from gravitational lensing (e.g.|Rusin et al"1l2Q03b 
iRusin & Ko chanek 2005; Jiang & K ochanek 2007|). In other 
words, the increase in the central density of the system ow- 
ing to the different degrees of dissipation in the two mergers, 
while apparent in the nuclear profile in Figure|2j does not rep- 
resent a sufficiently dramatic change in the total mass profile 
shape to significantly alter the virial constant (i.e. to signifi- 
cantly change the ratio M tme (r)/M^ yn (r)). That the estimator 
Md yn is similar across these extreme cases demonstrates that, 
despite the subtle difference in the strength of the central (dis- 
sipational or starburst) component in the stellar light profile, 
the homology assumption is valid in a purely empirical sense. 

The role of dissipation in driving tilt in the FP is not, then, 
to introduce substantial structural non-homology (in the sense 
of changing the profile sufficiently to alter the coefficient of 
the virial scalings) or to produce kinematic non-homology, 
but primarily to decrease the effective radius of the stellar 
light or mass distribution. At a smaller radius, the more 
stellar-dominated component of the system is sampled, and so 
when the traditional dynamical mass estimator is constructed, 
(Mdyn = ka 2 R e /G) and compared to the stellar mass, the ra- 
tio Mdvn/ M» is lower for systems wit h more dissipa t ion, a s 
argued in lRobertson et all d2006bl) and iDekel & Cox! d2006l) . 
Note that here a is either the central velocity dispersion or the 
average dispersion within R e (it makes little difference which 
we consider as long as we do so uniformly, since the disper- 
sion changes weakly with radius). 

Of course, / s b = 0.31 is a rather extreme case, but Figure[5] 
shows how the effective radius scales with gas content in an- 
other otherwise identical set of simulations of Milky Way- 
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FIG. 5. — Final stellar mass, dynamical mass Md yn = kR e a 2 /G, effective 
radius, and central velocity dispersion as a function of the dissipational or 
starburst mass fraction, for three sets of otherwise identical ~ L» disk-disk 
major mergers with varying initial gas content. Different line styles denote 
varying orbital parameters. The stellar mass is only weakly affected, but 
there is a strong scaling of size and dynamical mass (towards smaller sizes 
and correspondingly lower dynamical masses) with increasing dissipation (in 
the sense seen in Figure [3). 

like spiral mergers that span a range in dissipational frac- 
tion. The effective radii decrease systematically with dissi- 
pational/starburst fraction in a continuous manner. The ve- 
locity dispersions do become slightly larger as larger central 
mass concentrations are assembled, but the effect is much 
weaker (Alog a ~ 0.15dex, compared to A\ogR e ~ 0.8 dex, 
for / s b = 0-0.5). The result is that the dynamical mass esti- 
mated at R e decreases substantially with gas fraction - again 
owing to the effective radius of the stellar light being drawn 
further in - by AlogM^yn ~ 0.4 dex from low to high gas 
fraction. Meanwhile, the stellar mass is almost completely 
unchanged, since in any case most of the gas is eventually 
turned into stars, whether in the pre-merger disks or in the 
final starburst. So, driven by the increasing dissipation yield- 
ing smaller effective radii, the ratio Md yn /M* decreases by a 
factor > 2-3 as the gas fraction increases. 

6. ASSESSING THE ROLE OF DISSIPATION IN OBSERVED 
SYSTEMS 

Having illustrated how dissipation can alter the sizes of stel- 
lar spheroids and the ratio of dynamical to stellar mass in sim- 
ulations, we now proceed to apply our analysis to observed 
ellipticals. 

6.1. How Do Observed Sizes Scale with Dissipation? 

Figure [6] shows how effective radius scales with extra light 
mass /extra at fixed mass. We consider three mass bins, below, 
at, and above ~ M*. For each, we plot the residual R e relative 
to that expected for the given stellar mass, as a function of the 
fitted extra light frac tion. Specificall y, we determine (R e {M*)} 
from the sample of lShen et alJ (120031) . and take the ratio of the 
half mass radius of each system (determined directly from the 
light profile, or from the fits, it does not change the compari- 
son) to that value. Our mass bins are small enough, however, 
that this makes little difference compared to just e.g. consid- 
ering R e in a given bin. We show results separately for central 
cusp and core ellipticals, as well as our simulations, although 
the three distributions are similar. 

There is a significant trend: at a given stellar mass, systems 
with fractionally more extra light have systematically smaller 



R e {a also rises slightly, but the effect is weaker and there is 
more scatter at a given /extra. as expected owing to the role of 
the extended stellar and dark matter distribution in setting the 
central potential and a). In each case, the simulations and ob- 
served systems occupy a similar locus. We can also construct 
this plot with the starburst mass fraction / srJ of the best-fitting 
simulation as the independent variable (instead of the fitted 
extra light mass fractions /extra), an d find a correlation of the 
same nature. Given two progenitors of known size and mass, 
it is straightforward to predict the size of the remnant of a 
dissipation less merger, simply assuming energy conservation 
(see e.g. [B arnes 1988); in the case of a dissipative merger, 
ICovington et alJ (120081) use the impulse approximation to es- 
timate the energy loss in the gaseous component, followed 
by collapse in a self-gravitating starburst. This yields a de- 
tailed approximation as a function of e.g. initial structural and 
orbital parameters, but if we assume typical initial disk struc- 
tural scalings and parabolic orbits, it reduces to the remark- 
ably simple approximation 

7? e (dissipationless) 

e ™ i+CW/o) ' () 

where /o 0.25-0.30 and /? e (dissipationless) is the radius 
expected for a gas-free ( f s h = 0) remnant. (A similar estimate 
is obtained by Ciotti et al. (2007) assuming, instead of the im- 
pact approximation, a constant fractional dissipational energy 
loss). We plot this in Figure[6] with the scatter seen in the sim- 
ulations. At all masses, in both simulated and observed cusp 
and core ellipticals, more dissipational ellipticals are smaller 
in the manner predicted. 

This result directly implies that a structural difference - 
albeit not traditional structural non-homology - plays a key 
role in establishing the fundamental plane tilt. At fixed mass, 
smaller systems are so because a larger fraction of their mass 
is formed in a central dissipational starburst (to the extent that 
our extra light fractions recover this dissipational component). 
This dissipational starburst is compact, so even though the 
pre-existing stars are distributed to large radii, the effective ra- 
dius is smaller. In Paper II, we considered the light profiles of 
observed systems along the R e -f exaa correlation, at fixed stel- 
lar mass: it is clear that observed systems at fixed mass with 
the largest R e show profiles close to a pure Sersic law, with 
little evidence for any extra light component at > 0.01 R e (in- 
deed, they have /b < 0.02). Observed systems in this regime 
can still be classified as "cusp" ellipticals, but the cusps tend 
to be prominent at small radii and (in several cases) somewhat 
shallow, and contribute negligibly to the stellar mass. How- 
ever, moving to smaller R e and higher / ex tra at fixed stellar 
mass, deviations from a Sersic law at r <C R e become more 
noticeable. This is not to say that these deviations are univer- 
sal (that the extra light always takes the same shape/form), but 
there are increasingly prominent central light concentrations. 
If the systems were perfectly homologous (in the strict sense 
of the definition), there should be no differences, and there 
should be no trend whatsoever in Figure [6] 

Figure [6] also shows that there is little difference, in this 
sense, between the size scalings with dissipation of cusp and 
core ellipticals. We demonstrate this in Paper III, in both sim- 
ulations and observations. In short, this is expected: even if 
the systems with cores have expanded via re-mergers, they 
should (so long as there is not a wide range in number of re- 
mergers or strong systematic dependence of the number of re- 
mergers on starburst fraction at fixed mass) grow uniformly, 
and preserve these trends. We might expect some normal- 
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FIG. 6. — Effective radius R e relative to the median value for all ellipticals of the same stellar mass, as a function of our fitted extra light fractions (the 
empirical tracer of the dissipational/starburst mass fraction). We show simulated gas-rich merger remnants (top), observed cusp ellipticals and gas-rich merger 
remnants (middle), and observed core ellipticals (bottom); we use this point notation (key) throughout. We show this in three bins of stellar mass (relative to 
M» Ri 10" Mp, or M* = -21) Solid (dashed) lines show the mean (± 1 cr) correlation, following the analytic solution for dissipational mergers and fits to our 
simulation in Covington et al. 12008). We show the Spearman rank correlation coefficient and probability of the null hypothesis P m u (no correlation) in each 
panel. Simulations and observations exhibit the same behavior: systems with smaller R e at fixed mass have systematically larger extra light fractions (in the 
sense predicted by Figure|3]for the starburst/dissipational mass fractions). This implies that, at fixed mass, systems are driven along the fundamental plane by the 
relative amount of dissipation involved in their formation. This behavior is true regardless of cusp/core status. 



ization offset: if the mean size-mass relation after a gas-rich 
merger is a power law R e oc M", and two such systems with 
mass ratio / (where / is the mass ratio of the secondary to the 
primary) are involved in a dry merger from a parabolic orbit 
(and preserve profile shape), then the remnant will increase in 
radius by a factor (1 +/) 2 /(l +f 2 ~ a ) relative to the primary. 
However, it has also grown in mass, so compared to ellipticals 
of the same (final) mass, its relative increase in size is only 
(1 + f) 2 ~ "/(l + f 2 ~") . For observationally suggested values 
a = 0.56 dShen et al.ll2003l) . this predicts a - 30% (~ 0.1 dex) 
relative size increase for major re-mergers with mass ratios of 
1:3 through 1:1. This is easily dwarfed by the effects of dis- 
sipation seen in Figure[6] which can alter the sizes of systems 
by nearly an order of magnitude at fixed mass. That there 
is not a substantial normalization offset between the trends 
for cusp and core ellipticals in Figure [6] therefore implies that 
typical core ellipticals have b een involved in on ly a modest 
number of re-mergers (see also lCiotti et a l. 2007), but it is not 
a strong constraint. The important point is that observations 
demonstrate the dominant factor controlling the sizes of a el- 
lipticals is the amount of dissipation, even for systems which 



have undergone re-mergers. 



6.2. 



How Does This Change Spheroid Dynamical Masses? 

Figure [7] plots the quantity Md yn /M* as a function of the 
amount of dissipation, quantified by / ex tra- The expected trend 
bas ed on our fiducial cases in § [5] and the dissipational mod- 



els dOfiorbe et al.ll2005tlRobertson et al.ll2006bl:lDekel & Cox! 
2006) is borne out in our library of simulations over a wide dy- 
namic range with fairly small scatter, despite a wide range of 
total masses, orbital parameters, initial gas fractions, assump- 
tions regarding the gas pressurization, and other progenitor 
properties. The same trend obtains for the observed systems, 
whether cusp or core ellipticals. The scatter is significantly 
larger among the observed systems, but this is not surprising, 
given both the measurement errors involved and the possibil- 
ity for more complex growth histories that our simulations 
do not completely model. In a formal sense, the inverse cor- 
relation (i.e. smaller Mjyn/M* at larger dissipational fraction 
/extra) is highly significant: we plot the Spearman rank cor- 
relation coefficient p of each sub-sample, together with the 
significance of its deviation from zero. This P nu n represents 
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FIG. 7. — Ratio of dynamical to stellar mass within R e , as a function of extra 
light fraction derived from surface brightness profiles. Results are shown 
for simulated gas-rich merger remnants (top), cusp ellipticals and gas-rich 
merger remnants (middle), and core ellipticals (bottom). Dotted line shows 
an approximation to the mean simulation trend. We show the Spearman rank 
correlation coefficient (p) for each sub-sample and the probability P nu n of 
the null hypothesis that there is no correlation between Md yn /M* and /extra- 
Systems with larger degrees of dissipation (larger /extra) have larger ratios 
A^dyn/M* , with high significance, in the sense expected from Figure[3] 

the probability of the null hypothesis of no correlation be- 
tween Mdyn/M, and / ex tra, which is what would be expected 
if systems were perfectly homologous or if extra light did not 
contract systems and yield smaller R e . 

It is possible, in principle, that the observed correlation 
could be indirect. For example, if lower mass systems hap- 
pened to have higher /extra, but it had no causal connection to 
their lower Md yn /M». In Figure [8] we therefore consider this 
same comparison in narrow bins of stellar mass. We note that, 
while there is still some width to our bins, they are sufficiently 
narrow that we obtain the same answer if we assume that all 
systems in the bin have the same stellar mass. We also check 
a series of Monte Carlo experiments where we assume that 
Md yn /M* and / ex t m are both pure functions of Md yn or M*, 
but have no dependence on each other, and we find that the 
mass bins plotted in Figure [8] are sufficiently narrow that if 
this null hypothesis were true, then the maximal indirect cor- 
relation would have a significance of only p ~ -0.1 to —0.3 
(with P nu n w 0.2-0.8); so a significance P na a «C 0.1 is indeed 
meaningful. We caution that at a narrow range in stellar mass, 
there is naturally less observed dynamic range in dissipation 
fraction (as expected, since e.g. disks of a given mass have 
a reasonably narrow range of observed gas fractions), and the 
number of objects in each bin is smaller, so the significance of 
the inverse correlation in any one bin will be correspondingly 
reduced. However, we still see an inverse correlation in all 
bins, with significance stronger than expected from the null 



hypothesis. Near ~ M* where the sample size is large, the in- 
verse correlation is still highly significant. Further, while the 
significance in any one of the plotted bins may be lower ow- 
ing to the binning, the cumulative significance (i.e. likelihood 
of obtaining such consistent inverse correlations, considering 
each mass bin as an independent subsample) is actually quite 
high CP null « 10" 6 ). 

Figuresl9l&[T0lrepeat this test, using the starburst mass frac- 
tion / s b estimated from fitting simulations to observed light 
profiles as the proxy for the dissipational fraction. We obtain 
similar answers to our previous comparison using the fitted 
/extra, suggesting that our comparison is not peculiar to the ex- 
act estimator used, so long as it robustly recovers the physical 
dissipational component of the galaxy. Again, the cumulative 
significance of the inverse correlation - namely that Md yn /M* 
is smaller with increasing /b at fixed mass - is very high, 

Pnull~10- 8 . 

7. IS DISSIPATION NECESSARY TO EXPLAIN THE FP? 

These comparisons, while demonstrating that the dissipa- 
tional mass fraction does indeed correlate with the sizes and 
ratio Mdyn/M„ in observed ellipticals (and therefore that it 
could be the source of the tilt in the FP) do not necessarily in- 
dicate how much of the FP tilt derives from dissipation. To test 
this, we construct the fundamental plane in FigureQT]- specif- 
ically the correlation between Md yn and M* - in bins of similar 
fitted "extra light" mass fractions / extr a (i-e- we consider the 
FP for systems only with similar amounts of dissipation). Be- 
cause we are binning by extra light mass and still attempting 
to fit a correlation, we include both cusp and core ellipticals 
in the observed sample. Separating the two gives identical 
results, but the significance is reduced owing to the limited 
dynamic range from further splitting the sample at fixed / ex tra- 
We fit a power-law to the data (and separately, to our simu- 
lations) in each bin, of the form Md yn oc M\ +a (Equation [T), 
where a is the FP "tilt." Because neither Md yn nor M* should 
properly be considered an independent variable, we quote the 
results from fitting the least-squares bisector to the correla- 
tion. However, this makes little difference (the least-squares 
best-fit Md yn (M„) and M*(Md yn ) relations are both consistent 
and agree). We show this using the typical estimator Md yn 
here, but have also considered the true total enclosed mass 
Mtrue within R e and find identical results (as expected, since 
we show in §|5]and §|9]thatMd y n is a good proxy forM trae ). 

The results are unambiguous in both our simulations and 
the observed systems. The tilt in any given bin is negligible - 
within 1 a of a = in every case. In other words, at fixed / extra , 
we recover the virial correlation with constant Md yn / M» . Note 
that in several of these bins the dynamic range in mass is still 
large (^2-3 orders of magnitude in M„), and the number of 
observed systems is sufficiently large that if a substantial tilt 
like that observed for the entire population (a ~ 0.2) were 
present in the bin, we should see it. This is reflected in the 
quoted errors, which show that the tilt in most bins is incon- 
sistent with that observed for the global population at ~ 3 a 
significance (in the sense that it is much smaller, consistent 
with a pure virial correlation). If we combine the data (i.e. 
normalize out the mean M& ya /M* at each bin in / ex t ra and fit, 
or take the cumulative significance of the bins shown) we ob- 
tain a = 0.028 ±0.040, consistent with zero and ~ 4- 5 a be- 
low the observed tilt for the entire population. 

Figure Q~2] repeats this exercise, using the starburst mass 
fraction /b (estimated from the physical starburst mass frac- 
tions in the best-fitting simulations) as the proxy for dissipa- 
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FIG. 8. — As Figure \5\ but in narrow bins of stellar mass (range in M* for each bin shown in the upper right of each panel). Cusp (bottom) and core (top) 
ellipticals are shown. The same trend - towards lower M,j yn /M, is seen at fixed M t . The significance is lower in any one bin owing to the smaller number of 
points, but the cumulative significance is comparable to or higher than that in Figure[7] 
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FIG. 9. — As Figure[7] but using the starburst mass fraction (the physical 
mass fraction in our simulations, or that inferred from directly fitting simula- 
tions to the observed light profiles for the observations) instead of the purely 
empirical /extra as our estimator of the dissipational mass fraction. The results 
are similar. 

tional mass fraction, instead of the fitted extra light fraction 
/extra- In either case we reach identical conclusions for both 
the simulations and observations. The inferred tilt at fixed 
dissipational fraction is always consistent with zero, in some 
cases scattering to a < (i.e. the opposite sense of what is 
observed; but again these are all consistent with a = 0). In 



short, we have demonstrated empirically that without invok- 
ing some systematic dependence of dissipational fraction on 
mass, there is no FP tilt. Dissipation is observationally neces- 
sary to explain the FP tilt. 

Figure [13] also shows this for the effective radius projec- 
tion of the fundamental plane. We obtain identical answers in 
this case (for both simulations and observations) but note that 
the sense of the correlation is such that, if the baseline in R e 
is comparable to the scatter in observed points, severe biases 
can be introduced. We see this reflected in the fact that, in this 
case, fitting R e (R P Ted) versus R pTe ii(R e ) (where R pK d oc a 2 I~ l is 
the virial expectation) yield rather different slopes. Therefore, 
while we show our simulations (which have small dispersion 
in R e across sightlines), we refrain from fitting the observa- 
tions in such narrow bins in this particular projection of the 
FP. 

8. PREDICTING THE FP: IS DISSIPATION SUFFICIENT? 

8.1. Systematic Dependence of Dissipation on Mass: What Is 
Expected and Observed 

We have shown that Md yn /M* depends on the degree of dis- 
sipation, reflected in / ex tra or / s b, and that at fixed / ex tra or / s b 
there is essentially no FP tilt. This already implies that dis- 
sipation must be responsible for the majority of the tilt. But 
if we knew how much dissipation was expected as a function 
of galaxy mass (i.e. the mean expected mass fraction in a dis- 
sipational starburst, for mergers of systems of a given mass), 
we could entirely predict the FP and its tilt. 

Figure [14] shows the mean dissipational mass fractions of 
both cusp and core ellipticals, as a function of mass, derived 
in Paper II and Paper III. There is a clear systematic trend 
(discussed in those papers): lower-mass systems have sys- 
tematically higher dissipational mass fractions. To the extent 
that M<j vn /M* depends on dissipational fraction (/extra or /b), 
then, the existence of a systematic dependence of / ex tra and 
/b on mass will yield a systematic dependence of Md yn /M* 
on mass. We can use the observed dependence of the dis- 
sipational mass fraction on mass to predict and empirically 
estimate the amount of tilt in the FP contributed by system- 
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As Figure [8] but again using the physical starburst mass fraction instead of the fitted extra light mass fraction /extra- The results are similar 
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FIG. 1 1 . — The stellar-mass FP (i.e. M dyn — M* correlation), for simulated 
merger remnants (black solid) and observed ellipticals (open colored points), 
restricted to systems with a narrow range of extra light/dissipational mass 
fractions /extra (as labeled). In each bin of /extra, we fit the simulations and 
observations separately to a correlation of the form Md yn oc M' +a , where a 
quantifies the FP "tilt." The resulting tilts are shown (red upper and black 
lower values are for observations and simulations, respectively). The tilt is 
consistent with zero in every bin, and is lower than the expected tilt of the 
entire population (a ~ 0.2) by ~ 3cr - i.e. there is no "tilt" in either observed 
or simulated populations at the same, fixed value of /extra- 

atic trends in dissipation as a function of mass. For conve- 
nience, we noted in Paper II that the dependence of dissipa- 
tional fraction / S b on mass in cusp ellipticals (or observed gas- 
rich merger remnants) can be well-fitted by 



(./starburst) 



\M Q J 



(10) 



with (M , a) = (10 y ' 2±a2 M Q , 0.43 ±0.04) (shown in Fig- 
ure [141 . and a roughly constant factor ~ 1.5-2 scatter 
at each M*. In Paper III we demonstrated that an essen- 
tially identical scaling (a statistically equivalent (Mo, a) = 
(10 8 8±0 - 3 M Q , 0.35 ±0.05)) applies to core ellipticals, as ex- 
pected if cusp ellipticals are indeed their progenitors (since 
identical re-mergers will conserve dissipational or "extra 
light" mass fractions). 
We also demonstrated in Paper II and Paper III that this 
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FIG. 12. — As Figure fill but a gain using the physical starburst mass frac- 
tion / s (, instead of the fitted extra light mass fraction /extra- The results are 
similar in either case - there is a highly significant lack of "tilt" at fixed dis- 
sipational fraction, despite the observations and simulations still spanning a 
large baseline (~ 2dex) in mass. 
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FIG. 13. — The FP in narrow bins of / s t,, for our simulations, as Figure [121 
but showing an alternative projection of the FP (R e oc (cr 2 /j')' +a ). The result 
(no tilt at fixed f s ^) is the same; but in such a projection (with a small sample) 
small errors in R e will tend to amplify the observed tilt (and the physical 
meaning of the tilt is less clear). 
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FIG. 14. — Inferred gas content (dissipational/starburst fraction) of elliptical-producing mergers as a function of stellar mass. Initial gas fraction (top) and 
phy sical final starburst ma ss fraction (bottom) corresponding to the best-fit simulations to each observed system in the samples of Lauer et al. 12007a) (circles) 
and Kormendy et al. (2008) (stars) are shown, with the typical 25 — 75% allowed range (error bar). We s how results separately for cusp ellipticals (left), core 
ellipticals (center), and both together (right). Dashed (solid) line shows the fit to the data (Equation 1 10) in cu sp (core) ellipticals. Colored points with error 
bars indicate the mean (and rfc 1 a range in) disk gas fractio ns at the sam e stellar mass, at z = IBell & de Jo ng 200 1 ; Kannappan 2004 ; McGaugh 2003, blue 
diamonds, squares, and circles, respectively), z = 1 iShapley et al. 2005, green squares), and z = 2 )Erb et al. 2006:, orange circles). There is a clear trend of 
increasing dissipation required to explain elliptical profiles at lower masses (significant at > 8 cr in either core or cusp subsamples separately), in good agreement 
with the observed trend in progenitor disk gas fractions over the red shift range where cusp ellipticals are f ormed, and with what is invoked to explain the observed 
densities and fundamental plane correlations of ellipticals (e.g. Kormendy 1989; Hernquist et al. 1993). The best-fit trends in cusp and core populations are 
statistically identical: i.e. the dissipational/extra light component is preserved regardless of possible re-mergers. 



empirically measured systematic dependence of dissipational 
mass fraction on stellar mass agrees well with the observed 
dependence of disk galaxy gas fractions on mass. This is 
exactly what is expected if ellipticals are (at least origi- 
nally) formed in gas-rich, disk-disk mergers (regardless of 
whether or not they subsequently experience re-mergers). We 
show this in Figure [14] with the obs erved disk gas f r action s 
as a function of stellar mass from iBell & de Jon g (2001); 
iKannappan] d2004l); fMcGaughl d2005l) at z = 0, as w ell as at 
z. = 1 dShaplev et al.ll2005l) and z = 2 dErb et alj|2006h . At all 
these redshifts, the trend of gas fraction is similar, and the ob- 
served gas fractions evolve relatively weakly with redshift (by 
a factor ~ 1-2), bracketing the range of dissipational mass 
fractions observed in ellipticals at each mass. 

In other words, regardless of the exact times of formation, 
a systematic trend similar to that observed in the dissipational 
fractions of ellipticals is inevitable if disks are the progenitors 
of ellipticals. Note that even if there is a systematic depen- 
dence of the time of first gas-rich merger on stellar mass in 
this range, the mixing of formation redshifts would be less 
important than the mean dependence of disk gas fraction on 
mass (if we assume the stellar population ages date the merger 
times - which is an upper limit to the magnitude of this effect, 
since the stellar ages at least in part reflect the similar trend in 
stellar ages of disks as a function of mass - then the system- 
atic effects added by this age scaling are not important). The 
dependence of disk gas fractions on mass is therefore the pri- 



mary driver of the observed dependence of dissipational mass 
fraction on elliptical mass. 

Among other things, this implies that we could predict the 
dissipational mass fractions of ellipticals as a function of their 
mass, based on the observed disk gas fractions at each mass. 
We therefore pursue the following exercise: we consider the 
consequences of adopting the observed scaling of disk gas 
fractions with mass as our expectation for the dissipational 
fractions in ellipticals as a function of mass. Since we have 
just shown they are equivalent, we obtain identical results if 
we select the empirically fitted dissipational fractions of el- 
lipticals as a function of mass; in either case, the important 
thing is that there is a well-defined, physically motivated and 
observationally confirmed systematic scaling of dissipational 
fraction with mass. 

8.2. The Consequence for the FP 

Given the known scaling of disk gas fractions as a function 
of mass, we can then focus on the subsample of simulations 
that obey this scaling. That is, at a particular mass scale, we 
know the range of disk gas fractions seen empirically. We 
therefore consider only simulations within that range of gas 
fractions, as opposed to others in our library that are more gas 
rich or gas poor than typical observed disks of the same stel- 
lar mass. We have considered the effects of redshift evolution 
(using the mean age of ellipticals as a function of their stellar 
mass as a proxy for their merger time, we adopt the gas frac- 
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tion of disks of the appropriate mass at that redshift), and find 
it makes relatively little difference (although some secondary 
correlations along the FP relating to stellar populations and 
the age of systems will be discussed in § 111. 2\ . This restric- 
tion imposes a mean dependence of extra light or starburst 
mass fraction on stellar mass similar to that observed in ellip- 
ticals and anticipated from disk gas fractions. Given this con- 
straint, we can consider the FP and other correlations obeyed 
by just these objects: i.e. we can ask the question: What are 
the scalings obeyed by systems that have the same systematic 
dependence of dissipational fraction on mass as expected and 
observed in the real Universe? 

Figure Q3] shows the results of this exercise. Now that we 
account for the cosmologically expected and observed depen- 
dence of dissipational content on mass, the FP is "tilted." 
Fitting a power law, we obtain a value for this tilt of a = 
0.19 ±0.02. We compare this to the observed systems, plot- 
ting the data on the same footing, since we have constructed a 
cosmologically representative comparison sample in our sim- 
ulations. We obtain an observed tilt of a = 0.21 ± 0.05, the 
canonical value in the literature. More important (since any 
formal fitted tilt will vary if fit over different dynamic range), 
we compare the predicted scaling directly to actual ellipticals, 
and find that the simulations and observed systems occupy a 
statisti cally identical locus in M d vn - M* space. As a che ck, we 
fit the iKormendv et alj d2008) and lLauer et alj d2007al) sub- 
samples independently, and obtain consistent values for the 
tilt. 

We show the tilt again in terms of R e , and obtain a simi- 
lar result (for this projection, the formal tilt value is expected 
to be similar to that forMd yn -M*). Fitting the observed sam- 
ples, we find that they trace an identical locus - but the power- 
law fits in this case are more sensitive to dy namic range (ap- 
parent by t he fact that t he form al tilt from Kor mendv et alj 
(2008) and lLauer et all d2007al) are slightly discrepant, de- 
spite the points tracing each other where they overlap). As 
discussed above, caution should be used fitting in this FP pro- 
jection (here, the dynamic range is much larger than the fit 
errors in R e , so there are not large biases introduced in such 
observational fits, but in a more limited subsample, results fit- 
ted in this fashion are less robust than those fitted in the more 
physical Mdyn - space, where the axes are independent). 
Given the larger uncertainties in this case, the agreement is 
reasonable. Note that the absolute values of M* and Md yn and 
zero-points here are sensitive to our adopted virial coefficients 
and IMF assumptions, but the tilt and our meaningful compar- 
isons are not. 

Figure [16] compares this result with the samples divided 
into cusp and core ellipticals. For both classes of objects, the 
agreement with the simulations is good, and the two appear 
to trace a continuous FP correlation. This is not surprising: 
it has been shown that massive, boxy ellipticals with cores 
trace a continuous FP with less massive, disky, cuspy ellipti- 
cals (at least at the massive end, where most of the core popu- 
lation resides; see iGerhard et al.l 12001b Ivon der Linden et al.1 
120071) . and we have shown that the scaling of Md yn /M* 
with dissipational fraction (and systematic dependence of dis- 
sipational fraction on mass) are similar in both cusp and 
core ellipticals (and, in simulations in Paper III, we explic- 
itly demonstrate that this should be true for both the origi- 
nal gas-rich merger remnants and dry re-merger remnants) . 
Furthermore, a number of studies fcarjelato et al.1 119951: 
Dantas et ail 120031: iNipoti et al.l l200l iBovlan-Kolchin et all 
20051: iRobertson et al.l l2006bT) ~have shown that dissipation- 



less spheroid-spheroid re-mergers (popular as a mechanism 
for producing cores in ellipticals) tend to preserve the FP, pro- 
vided that the number of re-mergers is modest; we discuss 
these issues further in § [T0l 

Given that the FP tilt in Figure Q3] arises owing to dissipa- 
tion - where we have shown in Figure [TT] that at the same de- 
gree of dissipation, the virial correlation (no tilt) is recovered 

- it should be possible to explicitly factor out at least some 
of this tilt. Figure [17] attempts this exercise: we determine 
the best-fit bivariate function Md yn (M„, /extra)- For simplicity 
and to minimize the free parameters involved, we (for now) 
adopt a simple power-law parameterization Md yn oc M" fj? xm , 
and determine the best-fit parameters a and b such that the 
perpendicular scatter (i.e. scatter about the least-squares bi- 
sector) is minimized; we then repeat this using the estimated 
starburst mass fraction / s b instead of the fitted extra light frac- 
tion /extra- in both cases, and regardless of whether we fit to 
the observations or to the simulations, the best-fit correlation 
is essentially the virial correlation (i.e. awl) coupled with a 
dependence on / ex t ra or / sr) . In other words, the best-fit bivari- 
ate FP assigns the FP "tilt" entirely to the role of dissipation 

- at fixed degree of dissipation, there is no tilt. Explicitly in- 
cluding the degree of dissipation in this manner demonstrates 
that it is indeed "as fundamental" as any bivariate correla- 
tion (it is preferred with respect to a purely mass-dependent 
Miyn/M*). These fitted FP projections are nearly identical to 
the simplest predicted FP in dissipation-driven models (also 
shown in Figure [T7|>: namely, that Md yn oc M* F( f em:l ), where 
the function F(/ ex tra) oc Md yn /M» has the form seen in Figure[7] 
(note that this is more complex than a simple power-law). 

In principle, this explicit inclusion of an /extra or /b depen- 
dence should be able to account for some of the scatter in the 
FP, and therefore might provide a plane with smaller scatter. 
In practice, the scatter in e.g. Md yn (M*, / ex tra) is only slightly 
reduced relative to that in just Md yn (M„). The reasons for this 
are twofold. First, as discussed in § 18.41 the majority of the 
scatter in the FP owes not to different degrees of dissipation 
but to the combination of sightline-to-sightline variations, dis- 
persion in progenitor disk properties, and measurement errors. 
Second, as discussed in §[4] and Paper II, an estimator such as 
/b or /extra is of course not a perfect tracer of the true starburst 
mass fraction; with typical factor ~ 2 uncertainties. In prac- 
tice, then, most of the reduction in scatter from factoring out 
/extra is negated by scatter introduced in the variation between 
/extra and the "true" starburst mass fraction. The result is that 
the scatter in the FP in terms of M» and /extra is comparable to 
or slightly less than that in terms of M* alone. Owing to the 
uncertainties in estimating / ex tra, then, this FP representation 
is not necessarily a substantially improved observational pre- 
dictor of Md yn or e.g. R e at fixed a and M». However, the FP 
in these terms is certainly comparable to the "traditional" FP 
in how tightly it relates Md yn and M», and it is "more funda- 
mental" in terms of a physical explanation for the FP tilt. 

Because our primary focus in this paper is the understand- 
ing of the observed z = FP correlations, we reserve a more 
detailed modeling of the evolution of the physical FP (in terms 
of e.g. stellar and dynamical mass) with redshift (as e.g. typ- 
ical merger histories and progenit or gas content evolve) for 
future work (Hopkins et al. 2008e). For now, we simply note 
that the evolution in physical parameters considered here is 
predicted to be weak, in agreement with dire ct observational 
const raints from high redshift weak lensing (Hevma ns et alj 
2006) and optical studies (note that there is expected evolu- 
tion in optical bands owing to stellar population effects; these 
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FIG. 15. — The FP for all observed elli pticals (bottom), in the s tyle of Figu re l 1 1 \ (left) and Figure [T3l (rig/ir), with the measured tilt (black for all observed points; 
red and violet show the tilt fitted for the Kormendv et al. 12008) and Lauer et al. (2007a) subsamples, respectively). The standard tilt a ~ 0.2 is now apparent. 
We compare with the predicted correlation for our simulations (top), where we consider only simulations that match (within the observed factor ~ 2 scatter) the 
observed systematic correlati on b etween e.g. progenitor disk gas fractions and stellar mass, or (equivalently) between dissipational/starburst mass fraction and 
stellar mass (shown in Figure [T4V Accounting for the systematic dependence of dissipation on mass (owing to the dependence of gas content on mass in the 
progenitor disks), we obtain a FP tilt in good agreement with that observed, with comparable small scatter. 

ficients of a, R e , and I e are independently free in determining 
the best-fit projection of the observations. If, in principle, the 
FP were much more complex than a simple tilt (i.e. if M/L 
depended substantially on R e or a at fixed Md yn ), then the 
simulations might match the observations in Figures [T51I2TI 
but not in this representation. In fact, we find good agree- 
ment with the observations. This is not surprising, since it has 
been established that the FP can be re presented as a tilt un- 
der some homology assu mptions (e.g. Jorge nsen et al.|[l996l : 
iPadmanabhan et al.ll2004l) . 

Figure [19] plots projections of the FP; specifically, the 
size-mass relation and the velocity dispersion-mass (Faber- 
Jackson) relation. We show these correlations for the same 
simulations used to construct the FP in Figures [T31fT8l since 
we showed in § 16.11 that e.g. effective radii scale systemat- 
ically with dissipational fraction, a robust prediction of the 
size-mass relation requires that the simulations considered 
have appropriate dissipational fractions for systems of their 
stellar mass. We show sizes and velocity dispersions as a 
function of stellar mass because this is how these correlations 
are generally presented in the literature, but differences are 
small and the comparison with observations is similar if we 
consider them as a function of dynamical mass. 

We compare with the same correlations in the observed 
samples for which we consider the FP, and with the mean 
correlations measured for ellipticals in the SDSS. The sim- 
ulations agree well with the observations - both those that 
we consider in our FP analysis, and the trends in the general 
population (it is also reassuring that our observed FP samples 
obey the same projected scalings, with similar scatter). The 




FIG. 16 . — The observed FP, compared with our simulation predictions as 
in Figure [TSI but divided into core (top) and cusp (bottom) ellipticals. The 
simulation predictions agree well with the data in both cases, consistent with 
theoretical modeling that even a modest number of re-mergers do not move 
systems significantly off the FP correlation established by the original, gas- 
rich mergers. 

represent compli mentary constraints on elliptica l formation 
histories; see e.g.ldi Serego Alighieri et al J 120051: iTreu et al 



2005 



2007) 



Ivan der Wei et al. 2005; Ivan D okkum & van der M arel 
To the extent that these constraints agree with our ex- 
pectations based on empirical estimates for the evolution of 
disk gas fractions, they represent independent support for the 
scenario outlined here. 



8.3. Projections of the FP 

As anoth er check, we consi der our results in the observed 
FP space of iPahre et all d!998l) . in Figure [18] Here, the coef- 



Dissipation and the FP 



19 



12.5: AM = 0.07dex 

: AM(M„ «m;") = 0.08 dex 
12.0 : 

11.5 - 



5 

I 11.0- 

s 

° 10.5 - 
10.0 : 
9.5 



. ■ 



9.5 10.0 10.5 11.0 11.5 12.0 12.5 
0.99log(M,/lvy - 0.1 1log(f axlra ) + C 




12.5: AM = 0.17 dex 

: AM(M d «M!™) = 0.18 dex 
12.0 1 



11.5 r 
11.0 : 
10.5 : 
10.0: 
9.5 



: AM = 0.16 dex 



9.5 10.0 10.5 11.0 11.5 12.0 12.5 
O^logM/Ivy - 0.11log(f„ lra ) + C 



9.5 10.0 10.5 11.0 11.5 12.0 12.5 
logfM./Mg) + F(f eKlra ) 



12.5 : AM = 0.06 dex 

: AM(M a =m;*°) = 0.08 dex 
12.0 - 




12.5 


r AM = 0.16 dex 


= 0.18 dex - Vl 


12.0 


; AM(M dyn o=M!*») 


11.5 




e 


11.0 






10.5 






10.0 


'■ 




9.5 







9.5 10.0 10.5 11.0 11.5 12.0 
O.gSlogfM./Me) - 0.21 log(f J + C 



12.5 



9.5 10.0 10.5 11.0 11.5 12.0 
0.89log(M,/Me) - 0.41log(f sb ) + C 



12.5 




10.0 10.5 11.0 11.5 
log(M,/M 6 ) + F(fJ 



FIG. 17. — A "more fundamental" representation of the FP in simulations (left; points from Figure [T"5l and observations (center, right): Mjy,, as a bivariate 
function of M* and fitted extra light fraction (/extra, top) or inferred starburst fraction (/ s b, bottom). The best-fit correlation recovers M(j yn ex M„ with the residual 
dependence - i.e. the "tilt" - owing entirely to the degree of dissipation. Compared to the FP in terms of just M^ yn and M„ , the scatter is slightly reduced (but not 
much, owing to the uncertainties introduced in estimating /extra or /sb)- Adopting Mi yn oc M* with the additional functional dependence of Afd yn /M* on /extra 
from Figure |7]yields a similar correlation (right) with further slightly reduced scatter. 
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FIG. 18. — Observed FP, in the best-fit projection from Pahre et al. 1 1998); 
i.e. freeing the exponents in R e ono^ {] ■ Open points a re th e observed sys- 
tems, filled points are our simulation results from Figure [TH If the observed 
FP was more complex than a tilt (e.g. if M* depended substantially on R s or 
cr separately at fixed Md yn ), our results would disagree here. The agreement 
is good, as expected based on observations that the FP is, essentially, a tilted 
virial plane. 

scatter in these correlations is larger than that between Md yn 
and M*, since at fixed Md yn there is a tradeoff between R e and 
g. It is also interesting to note that the trends show signifi- 
cant curvature - their combination (Mdy n -M*) is reasonably 
approximated by a power-law scaling, but either the R e -M* 
or er-M* slope will vary with the fitted mass interval (for this 
reason, we refrain from quantifying a power-law fit to either 
the simulations or observations, but we note that if we do so 
over the same mass interval for both, we obtain the same re- 
sult). 

This has been discussed in the liter ature, for both the 
R e — and a — M* relations (see e.g. lLauer et all l2007bh 



FIG. 19. — Projections of the FP. We show the size-mass (top) and velocity 
dispersion-mass (Faber-Jackson; bottom) correlati ons, for the same subset of 
simulations used in our FP predictions in Figure 1151 (left) and observations 
(right). The observed mean relations for SDSS ellipticals (solid lines; with 
~ 95% ran ge as dashed lines) are shown over the r ange of data used in the ob- 
served fits (Shen et al. 2003; Gallazzi et al. 2006, for size-mass and velocity 
dispersion-mass, respectively). Note the curvature in both relations, owing to 
the dependence of dissipation on mass and its role in setting both R e and a. 
The simulations and observations agree with each other and with the mean 
correlations observed in much larger samples; i.e. once systems with the ap- 
propriate range of dissipational fraction for their mass are selected, they fol- 
low the observed size-mass and other correlations. The relations would not 
agree if we used systems with significantly different dissipational fractions as 
a function of mass. 



Ivon der Lin den et al. 2007t iDesroches et al.ll2007l) (note that 



even iShen et al.1 (120031) see tentative evidence for this effect 
in their lowest-mass bins, but their dynamic range at the low 
mass end is limited). We note here that such curvature can 
arise even from pure gas-rich merger remnants, owing to the 
dependence of the amount of dissipation on mass and cor- 
responding scalings of size and velocity dispersion. This is 
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supported by the fact that both our cusp and core sub-samples 
independently exhibit similar curvature. It has also been ob- 
served that non-BCG and BCG galaxy samples show similar 
curvature (despite the latter being much more likely to hav e 
experienced dry mergers; see e.g. Ide Lucia & Blaizoti 12007). 
This is not to say that re-mergers will not introduce curvature 
in these relations or move systems considera bly with respect 
to them (see e.g. iBovlan-Kolchin et al.l 120061) . but rather to 
emphasize that dry re-mergers are not the only source of cur- 
vature in the relations, and curvature should not necessarily 
be interpreted as evidence for dry mergers (further note that 
systems cannot move too far off their initial correlations, as 
they are constrained by the observed scatter). 

It is well-established that these correlations are different 
from those obeyed by spiral galaxies. If we do crudely 
approximate the correlations as power-laws (for illustrative 
purposes) over the dynamic range sho wn in Figure [T9l we 
obtain correlations R e oc M^' 56 (as in IShen et al] l2003h and 
a ocM^ 25 (see e.g. iBernardi et al.1 [2003). also yielding R e oc 
0-2.0-2. 2^ Compare these to the observed relations in disks: 
R„ oc M° 3 °-°- 3 \ V max oc Mg 27 " 030 , and R, oc V±t h2 (see e.g. 
Persic etalJ 1 1996bt ICourteau et ail 120071: lAvila-Reese et all 
2008). These differences relate to the tilt of the FP and the 
fundamental distinctions at issue here: at high masses, ellipti- 
cals and disks have similar sizes, but at low masses, ellipticals 
are much smaller. The difference, in our simulations, arises 
because of changes in the degree of dissipation, reflected in 
the extra light content. In principle then, it should be possi- 
ble to consider the resulting correlations between R e and M* 
or a in terms of the combination of their dependence on that 
variable and the additional dependence on / extra or / s t,. 

Figurel20lshows the results of this exercise: we consider R e 
as a bivariate function of combinations of M*, a, and /extra- 
The correlation R e (M*, a) represents the "traditional" FP cor- 
relation - that the correlation is not simply R e oc M»/<7 2 re- 
flects the "tilt" of the FP. If instead we fit R e to a combina- 
tion of M» and / ex tra> we obtain a fit of similar quality (i.e. 
this combination is an equivalent predictor of R e ), but with 
the dependence R e oc f^?. This is an equivalent FP 
in terms of its usefulness as a predictor of R e , with a scatter 
~ 0.15-0.20dex. More interesting, the explicit dependence 
of R e on M* fitted in this manner is identical to that observed 
in disks (R e oc M° 3 ), with the remaining dependence owing to 
the dependence of R e on /extra- Given the cosmological scal- 
ing of / ex tra with M*, we obtain the steeper R e oc Mj 6 relation 
for ellipticals when / ex t ra is ignored, but at fixed / ex tra or / s b 
- i.e. fixed degree of dissipation - a relation R & oc M ( J 3 , simi- 
lar to that observed i n disk galaxies, is obser ved (in line with 
the predictions from Robertson et al. 2006b). We obtain the 
same results for the size-velocity dispersion relation: a best- 
fit R e oc a h33 f^f (or R e oc a 0M f^f ); the dependence of 
/? e on a at fixed / ex t m or / S b is very similar to the dependence 
of R e on rotational velocity V max in disks (to lowest order, our 
simulations and other numerical experiments demonstrate that 
a traces the pre-merger V max , modulo a roughly constant nor- 
malization offset owing to e.g. the profile shape and kinemat- 
ics). In other words, the difference between the observed size- 
mass-velocity dispersion relations in disks and ellipticals can 
be entirely accounted for by dissipation. 

8.4. The Small Scatter in the FP 

A final requirement for any model of the FP is that it ac- 
count for the small observed scatter. Analyzing the relations 
predicted by the simulations, which broadly sample orbital 



parameters and stellar masses, and fully cover the observed 
range of dissipational fractions at each mass, we can see di- 
rectly in Figures [T51I2T1 that there is small scatter. Formally, 
we find a nearly symmetric la scatter of ~ 0.07- 0.08 dex in 
M dyn (M*) or ~ 0.065-0.075 dex in M*(M dyn ) (0.055 dex scat- 
ter perpendicular to the best-fit correlation - or equivalently 
in the R e projection of the FP - depending weakly on how 
we define our expected ffmAM*) and whether we use M tme 
orMdyn). This is the scatter obtained in these correlations us- 
ing the median values of M* and Md yn for each galaxy (across 
~ lOOsightlines): if we include the sightline-to-sightline vari- 
ance in the simulations, we obtain ~ 0.14 dex total scatter in 
Md yn (M*) (0.09 dex scatter perpendicular to the correlation). 
This is still less than that observed in each case: for the ob- 
served samples herein, we find ~ 0.18 dex scatter in Md yn (M*) 
and ~ 0.15 dex scatter in M*(Md yn ); i.e. a perpendicular scat- 
ter of 0.12 dex, with additional observational errors in R e and 
M* (estimated from our experiments in Paper II) likely con- 
tributing most of the difference (for typical ^0.1 dex stellar 
population model errors in M», this yields roughly the ob- 
served scatter). 

The origin of the small intrinsic scatter in the FP is straight- 
forward to understand in this scenario. Some intrinsic vari- 
ation will come from the scatter in the total baryon-to-dark- 
matter content of the progenitor galaxies. However, obser- 
vations of the baryonic Tully-Fisher relation suggest that the 
scatter in M* at fixed maximum circular velocity (a proxy 
for halo mass) is sm all, m 0.1 dex (at lea st over the mass 
range of interest here: lBell & de Jo ng 2001). Even consider- 
ing the total baryon-to-halo mass ratios, which extend well be- 
yond R e and are therefore of l ess interest, observations imply 
quite s mall scatter ~ 0.15 dex dWang et al.l2006t IConrov et al.l 
2006; Weinma nn et al.l l2006h . This effectively subsumes a 
number of quantities, including e.g. scatter in the initial radii 
of the disks at fixed mass, in halo concentrations, and in other 
parameters (note that the scatter in any one of these quanti- 
ties is not important, given that the final Tully-Fisher relation 
scatter is small). The contribution of this initial scatter will 
actually be reduced after a gas-rich merger, because the quan- 
tity of interest is the ratio [M»/2+M dm (< R e )]/(M*/2). Since 
the merger channels gas into the center and raises the cen- 
tral density, the final effective radius is smaller (explaining 
the smaller sizes of ellipticals relative to spirals at low mass). 
Thus, while M* is just the baryonic mass of the progenitor, 
the contraction of the remnant makes the relative value of 
Mdm(< Re) smaller, reducing the importance of initial scatter 
in Md m /Mbai y on to the final FP scatter. Analysis of the simu- 
lations suggests that a realistic initial ~ 0.1 dex scatter in the 
baryonic Tully-Fisher relation contributes < 0.04 -0.05 dex 
scatter directly to the final FP prediction. 

However, there will be additional intrinsic scatter con- 
tributed by different degrees of dissipation in systems of the 
same mass. Figures|6]&|8]shows that at fixed stellar mass, sys- 
tems with different degrees of dissipation have correspond- 
ingly different effective radii and ratios Md yn /M». However, 
the scatter in the expected (and observationally inferred) dis- 
sipational fraction at each mass is small, a factor < 2, and 
the dependence of Md yn /M» (see Figure |7]l on dissipational 
fraction over this relatively narrow range is not strong. Com- 
bining this scatter with the dependence of Md yn /M*, we ex- 
pect it to contribute about ~ 0.06 -0.08 dex of intrinsic scat- 
ter. Together, this yields the relatively small ~ 0.08 dex in- 
trinsic scatter in Md yn (M*) (before sightline-to-sightline and 
measurement errors). 
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FIG. 20. — FP predictors of R e . Left: R e as a function of the "traditional" global parameters M„ (stellar mass) and <r (velocity dispersion). Simulations (the 
same as used for the predicted FP in Figure [TSl and observations are shown (as in Figure[7}. The scatter in R e in this projection is ~ 0.2dex. Center: R e as a 
bivariate function of M, and fitted extra light fraction (/extra, top) or inferred starburst fraction (/ s b, bottom). This is an equivalent or better observational predictor 
of R e (scatter ~ 0. 15dex). Right: R e as a function of <r and /extra (top) or / s b ( bottom). Again, this is an equivalent predictor (scatter ~ 0.2 dex). Ignoring the 
effects of dissipation, we recover the R e -M* and a-M t relations in Figure [T9l with very different slopes from the relations for disks (comparison shown here). 
However, with the effects of dissipation explicitly included here, the relation between R e and Af» or R e and tr is equivalent to that obeyed by spirals - i.e. the 
difference between the "projected" elliptical and spiral scaling relations is entirely attributable to the degree of dissipation (/extra or / s b). 



9. TOTAL VERSUS DYNAMICAL MASS: INVARIANCE OF THE FP 

In Figure [21] we repeat our experiment from Figure Q3J 
i.e. we construct the FP from the simulations, given the ob- 
served and expected dependence of dissipational fraction on 
mass. However, instead of plotting the traditional dynami- 
cal mass estimator Md yn oc a 2 R e , we plot the true total pro- 
jected mass (stellar plus dark matter plus gas, although gas 
is generally negligible) within R e , which we extract directly 
from the simulations and refer to as M tme . The predicted FP 
is nearly identical to that predicted using the Md yn estimator 
- i.e. the FP predicted does, in fact, arise from a change of 
the dark-to-luminous-matter ratio within R e , owing to dissi- 
pation contracting the stellar effective radius (i.e. a change 
in Mtme/M* with mass), rather than to traditional structural or 
kinematic non-homology (which would imply a constant ratio 
of Mtrue/M*, with a changing apparent mass from the dynam- 
ical mass estimator; i.e. a varying M true /Mdy n with mass). 

This is despite the systematic variation of extra light frac- 
tion with mass (reflecting the systematic dependence of dis- 
sipational or initial disk gas fraction). In other words, while 
there is technically a subtle non-homology implicit in the fact 
that the dissipational fraction depends on mass (and, as a re- 
sult, the fitted / extr a and / s h change with mass, whereas they 
would be constant for true perfectly self-similar systems), it 
does not contribute any significant or observationally mean- 
ingful structural or kinematic non-homology. The "homology 
assumption," namely that the ratio of the dynamical mass es- 
timator to the true enclosed mass within the effective radius 
(Mtrue/Mdyn) is constant, is predicted to hold, while the FP is 
satisfied. 



Figure [22] demonstrates this explicitly: we plot the ratio of 
true enclosed mass to the dynamical mass estimator, specif- 
ically for the same simulations we used to predict the FP 
and its tilt. The ratio M tme /Md yn is constant despite chang- 



ing dissipational fractions; fitting a power law M tlue oc M^ , 
we obtain 5 = -0.01 ±0.01 (both insignificant and completely 
negligible compared to the observed value of the FP tilt, 
a » 0.2); i.e. the "homology assumption" is valid in the ob- 
servational sense as applied to these simulations, and there is 
no significant implied systematic structural or kinematic non- 
homology (in the traditional sense of the term) as a function 
of mass. 

We compare these predictions with observed "true" en- 
closed masses within the effective radii of observed sys- 
tems, using two different approaches for estimating the true 
enclosed mass (with out explicitly in v oking the homology 
assumption). First, Cappellari et al. (2006) estimated to- 
tal masses from three-dimensional Schwarzchild modeling 
or two-dimensional Jeans modeling of local ellipticals with 
two-dimension al velocity field information from SAU RON 
maps (see also [E msellem et al. 20Q3 iMcDer mid et al] 120061 : 



ICappellari et al.l 120071) . and~used stellar population models 
to estimate stellar masses (we correct these for the choice 
of IMF, but otherwise do not modify them). We also have 
considered integ ral modeling masses from other, independent 
sources (e.g. Ivan der Marell Il99lt iKronawitter eTafl 120001: 
Haring & Rix 2004) and obtain identical results in each case. 

Figure[2T]plots the FP of these objects; i.e. the true enclosed 
mass determined in this manner as a function of stellar mass. 
The tilt, in agreement with the simulations, is nearly identi- 
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FIG. 21. — Predicted and observed FP (as in Figure \l5[ , but using the 
true total enclosed mass M llue (stellar plus dark matter ) within R e in s tead o f 
Md yn . We compare with the observed systems from Cappella ri et al.l 120061) . 
who use integral Jeans/Schwarzchild modeling to determine a "true" mass 
independent of a homology or constant M/L assumption. A nearly identical 
tilt is predicted and observed - i.e. the predicted and observed FP is really a 
reflection of the tilt in M lrll e/M» within R e , not a tilt in Mtms/M^m. 
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FIG. 22. — Top: The correlation between the dynamical mass estimator 
Mdyn = 3.8 a 2 R e /G and the tru e tot al mass Mtme (stellar plus dark matter) 
enclosed within R e (as in Figure |2"T\ . for the same simulations shown in our 
predicted FP (Figure[l5j. Dashed line is a fit of the form Mt rue oc with 
the best-fit value of S shown: the best fit is indistinguishable from the Mtrue = 
Md yn l ine (i.e. has no "tilt"). Botto m: Same, for the obs erved systems from 
Figure |2T1 using integral modeling I Cappellari et al. 2006, magenta squares) 
or lensing mass maps (Bolton et al. 2007, 2008, orange triangles) to estimate 
the true mass Mtr ue . The best-fit coefficients M truc oc are shown for each 
observed sample. The simulations we use to predict the FP, and the observed 
systems on the FP, trace a nearly identical M tmc oc <r 2 R e relation consistent 



cal to that given bythe dynamical mass estimator M dyn . This with the traditional observational definition of the "homology assumption." 

is clear in Figure [22] where we compare the predicted cor- 
relation between Md yn and M true with these observations (we 
apply the same definition of Md yn to the observed objects as 
to the simulations, nam ely Md yn = 3.8R e a 2 /G). As noted by 
Capp ellari et al.l (120061) . there is a tight correlation between 
Mdyn and without any significant deviation in this space 
from the homology assumption (5 = +0.012 ± 0.019). Not 
only does this tight proportionality (5 w 0) agree with our sim- 
ulations, but the normalization fitted to the observed sample, 
i.e. mean ratio M true w 3.8R e a 2 /G, agrees with that predicted 
by the simulations over the observed mass range to better than 
0.03 dex, further suggesting that the structural properties and 
profile shapes of the simulations are in good agreement with 

those obse rved. 

Second, Bolt on et~afl (120071 120081) have performed a sim- 
ilar exercise, but using instead the enclosed masses deter- 
mined from gravitational lensing. We compare their FP (i.e. 
M true (M»)) in Figure [2T| and again find similar tilt to that 
predicted by the simulations (using the true enclosed mass 
within R e as M tme ) and to that obtained using Md yn instead 
of M trae . Figure [22] shows that, for these objects as well, 
Mtme oc M dyn (8 = -0.004 ± 0.034 » 0), with again a normal- 
ization (mean M tl - ue /Md yn ) in agreement with the simulations 
to within 0.02 dex. Note that the authors restricted themselves 
to true masses and dynamical masses within R e /2; if we use 
their lens models to correct the true masses to R e and take 
the corresponding dynamical masses at that radius, we obtain 
identical results. 



The simulations reproduce the observed, nearly exact pro- 
portionality between M tlue and Md yn , and therefore show the 
same FP behavior as has been seen in the observations: 
namely that the FP tilt remains similar regardless of whether 
Mdy n or M trlle is considered. This demonstrates the point from 
§0 that non-homology (in the general, observationally mo- 
tivated sense of different M true /Md y n as a function of galaxy 
properties) is not a significant contributor to the FP tilt. 

10. RE-MERGERS AND FP EVOLUTION 

The predictions we seek to test, and the observed depen- 
dence of FP tilt on the degree of dissipation which we have 
demonstrated, hold in the same manner for both cusp and core 
elli pticals. However , a number of observed properties (see 
e.g. lFaber et al.l l997) suggest that core ellipticals may be the 
products of subsequent "dry" or spheroid-spheroid, relatively 
gas-free re-mergers of (generally lower-mass) cusp ellipticals 
(formed directly in gas-rich mergers). If true, do we expect 
the level of agreement we see? 

The answer is generally yes, provided that the number 
of re-mergers is no t large. Numerica l experiments (e.g . 
I Capelato et all 119951; [ Dantas et al. 120031: iNi poti et al.| 120031: 
Bovlan-Kolch in et al.l l2005: Rob ertson et al. 2006b) indicate 
that remnants of dissipationless mergers of systems which be- 
gin on the FP tend to remain on the FP. Even if the re-merger 
introduces no new tilt, its effects would be small. To lowest 
order, since the gas is mostly exhausted, the re-merger will be 
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dissipationless, and merging two identical systems M\ andM2 
on a parabolic orbit, one expects their profiles to be roughly 
preserved. This leads to the energy conservation equation 

E f = k(M 1 + M 2 ) a) = E t = kM\ a\+kM 2 a\ (11) 

where <r/ is the velocity dispersion of the final remnant, and 
k is a constant that depends on the shape of the profile. From 
this, one obtains the general rule that a major 1:1 merger will 
approximately do uble R e while preserving a, doubling both 
M dvn and M„ (e.g. lHausman & Ostrikerill978UHernquist et all 
fl993h . 

Relative to the FP scaling, Md yn oc M' 2 , the remnant is now 
~ 0.06 dex below the FP expectation - only ~ 0.3 -0.5 a (a 
being the observed FP scatter). So even if the system increases 
its mass by a total factor of ~ 4-5 via dry re-mergers (i.e. has 
~ two equal-mass re-mergers or ~ 5-6 more typical 1:3 mass 
ratio re-mergers), it will move by only ~ 1 a with respect to 
theFP. 

Compare this to the impact of dissipation, which as we 
demonstrate in § [6] can change the effective radius and ratio 
Mi yn /M* (at fixed stellar mass) by nearly an order of magni- 
tude. We therefore expect that, unless a system has experi- 
enced an extreme number of major dry re-mergers (expected 
only for the rarest, most massive central cluster galaxies), the 
degree of dissipation (i.e. properties we study herein) should 
be the dominant factor determining the effective radii, ratio of 
dynamical to stellar mass, and loca tion with respect to the FP . 

Furthermore, it has be en noted dBovlan-Kolchin et al.l2005l 
20061 ICiotti etai] l2007) that re-mergers can slightly increase 
MAy n /M*, essentially causing systems to move nearly paral- 
lel to the FP. We consider a sample of re-merger simulations, 
described in § [2] in order to demonstrate this in Figure [23] 
Specifically, these are identical (mass ratio 1:1) re-mergers 
of remnants of previous gas-rich, disk-disk mergers, merged 
in various orbital configurations. The details are discussed 
in § E] but we note that although we plot only identical re- 
mergers here (for illustrative purposes), the results scale ap- 
propriately for different mass ratios and various mixed en- 
counters (i.e. merging different spheroid remnants of similar 
mass). Plotting both the pre and post-remerger systems on 
the FP in Figure [23] we see that they move nearly parallel 
to the fitted relation (a « 0.2). Quantitatively, the median 
increase in dynamical to stellar mass ratio is ~ 0.06 dex, pre- 
cisely as needed to move the system along the FP (to offset the 
~ 0.06 dex offset from the FP estimated above, if Md yn /M* 
were exactly conserved in a re-merger). 

This effect is subtly related to e.g. changes in the orbital 
isotropy and kinematics of the remnant in a re-merger, but 
primarily owes to a real (albeit small) physical increase in the 
ratio of enclosed dark matter mass within the stellar effec- 
tive radius (the progenitors and remnants obey the "homol- 
ogy assumption" as in Figure [22] and we see an almost iden- 
tical effect plotting the true enclosed mass instead of Md yn ). 
Considering the initial and final distributions of stellar and 
dark matter particles as a function of their initial radius, it is 
straightforward to understand this effect. 

Despite the nearly uniform inflation of the light profile in 
a re-merger, there is substantial scattering of stars (and dark 
matter): i.e. although stars at some initial r, will be, at the end 
of the merger, at a median r/ « 2r,, the distribution of final 
radii will be approximately lognormal with scatter w 0.4 dex 
(see Paper III). The integrated stellar mass will be constant, 
but such a scattering will tend to broaden the dark matter 
distribution in a way that slightly increases the dark matter 
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FIG. 23. — Top: Impact of subsequent s pher oid-spheroid re-mergers on the 
FP. We show the observed FP from Figure [T"5l (black line), with the initial and 
final positions of si mulated systems before and after a major re-merger (point 
style as Figure [T51 lines connect each re-merger to its progenitor). Bottom: 
The cumulative distribution in relative dynamical to stellar mass ratio before 
and after a major re-merger, in our entire ensemble of re-merger simulations. 
Vertical dashed (dotted) line shows the median (Her) change in Md yn /M». 
Re-mergers tend to slightly increase Mj yn /M», by almost exactly the amount 
needed to move parallel to the FP (0.06 dex), but with non-negligible scatter 
(~ 0.065 dex). 

mass within the stellar R e - i.e. scatter some dark matter from 
near the effective radius of the halo to both smaller and larger 
radii (total dark matter mass is conserved, but the densities are 
slightly lower near the halo effective radius - which is much 
larger than the galaxy effective radius which concerns us here 
- and slightly higher at much smaller and larger radii). This 
will slightly raise the central dark matter density, in a way 
consistent with the observed (weak) change inMd yn /M* in re- 
mergers (effectively mixing slightly more dark matter into the 
stellar distribution). The net effect of this is that re-mergers 
move systems even less with respect to the FP, and dissipation 
is still the dominant factor setting Md y n/M*. 

Figure [23] also demonstrates that although re-mergers tend, 
in the mean, to increase Md yn /M* such that they move parallel 
to the FP, there is significant scatter in the change in M^n/M* 
in re-mergers (in detail, the median increase is 0.06 dex, but 
there is a 1 u scatter of ~ 0.065 dex about this median). We 
therefore expect that growing a population of ellipticals by a 
factor ~ 2 in re-mergers (equivalent to a single 1 : 1 or three 1:3 
re-merger for most of the objects in the population) will con- 
tribute ~ 0.06 dex intrinsic scatter in Md yn (M*). Compared 
to the intrinsic scatter (in terms of sightline-averaged quanti- 
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ties) of ~ 0.08 dex, this is significant, yielding a final intrinsic 
scatter of ~ 0.10 dex. However, after accounting for the ad- 
ditional ~ 0. 1 dex sightline-to-sightline scatter, and observa- 
tional effects, the contributed scatter is negligible compared 
to the final observed ~ 0.1 8 dex scatter in Md yn (M*). That 
the observed scatter is not much larger does suggest there has 
not been dramatic re-merging, but the constraint is weak - 
it would take ^4-5 major (mass ratio 1:1) dry mergers (or 
~ 10 more likely 1:3 dry mergers) to noticeably increase the 
observed FP scatter at the massive end. 

11. DISCUSSION 

We have developed and implemented a set of observa- 
tional tests of the theoretical proposal that different degrees 
of dissipation are respon sible for the tilt i n the f undamental 
plane, ( suggested in e . g.lDjorgovski et alJ ( 1988h : iKormendyl 
(119891) : iBender et al.1 j 19921): ICiotti et al l (Il996h and devel- 
oped i n numerical simul a tions inlBekkil d!998 ): Onorb eet al.l 
d2005l) : iRobertson etail d2006bl) and iDekel & Cox! d2006l) ). 
With measurements of the surface brightness profiles of el- 
lipticals of sufficient accuracy and high dynamic range, we 
demonstrated in Paper I-Paper III that it is possible to design 
an empirical decomposition which reliably separates the dis- 
sipational (original merger-driven starburst) and dissipation- 
less (scattered stars from the original stellar disks) compo- 
nents of galaxy light profiles. Applying this to a large sample 
of observed ellipticals, we study here how their FP correla- 
tions relate to the mass fractions in the dissipational (or "extra 
light'Vstarburst) component. 

We show that systems with larger dissipational fractions 
have smaller effective radii for their stellar masses and lower 
ratios of total (stellar plus dark matter) mass to stellar mass 
within their effective radii (equivalently, lower Md yn /M*) - 
i.e. that dissipation can a nd does move ob j ects on to the FP, in 
the manner predicted by Robertso n et al.l d2006bl) . More im- 
portant, we show that without the effects of a systematic de- 
pendence of the degree of dissipation on mass, there is no sig- 
nificant FP "tilt". In other words, systems with the same dis- 
sipational/starburst fractions, even over a wide range in mass 
(~ 2 - 3 dex), exhibit a constant ratio M true oc M„ . Dissipation 
not only can tilt systems on the FP, it is required to explain 
the FP tilt. 

The tilt originates because of an expected systematic de- 
pendence of dissipational content on mass, arising because 
lower-mass disks (the progenitors of lower-mass ellipticals in 
this scenario) tend to be more gas rich. Therefore, lower-mass 
systems will, on average, undergo more dissipational merg- 
ers. This appears to be true even out to redshifts z ~ 2 - 3, 
so that the exact formation times of the systems are not espe- 
cially important. Quantitatively, if we convolve the expected 
dependence of dissipational fraction on mass (from e.g. the 
dependence of gas fraction on disk mass, or from the clear 
mean systematic trend we see in the fits to the observed el- 
lipticals) with the dependence of Md yn /M» on dissipational 
fraction, we obtain a prediction for the FP which is "tilted" in 
a manner similar to that observed, Md yn oc M* 2 . 

Of course, merger-induced starbursts may not be the only 
source of dissipation; for example, stellar mass loss may re- 
plenish t he gas supply and lea d to new dissipational bursts 
(see e.g. lCiotti & Ost riker 2007). Moreover, the merger his- 
tory and series of induced dissipational events may be more 
complex than a sing le or c ouple of ide alized major mergers 
(see e.g. iKobavashil r2004t iNaab et alJfeOOTl) . For our pur- 
poses, however, all dissipational star formation will appear 



similar when observed and have the same effects - we are 
essentially measuring the integrated amount of dissipation, 
and testing the idea that the presence of this dissipation can 
explain the difference between spheroid and disk scaling re- 
lations. In any case, the result is the same, and the agree- 
ment with disk gas fractions over the redshift range of interest 
suggests star-forming galaxies are viable candidates for the 
ultimate progenitor systems of ellipticals (however complex 
the details of the morphological transformation may be). In 
short, the systematic dependence of the degree of dissipation 
on mass in elliptical formation is both necessary and sufficient 
to explain the FP tilt. 

11.1. Possible Sources of Hit: Why Dissipation is the Only 
Viable Possibility 

This conclusion is, to our knowledge, consistent with all 
existing observational constraints on the FP. We note that 
repeating our experiments with alternative IMFs affects 
only the normalizations of the quoted relations; consider- 
ing alternative stella r evoluti on mo d els su ch as those of 
iLeitherer et al.1 d 19991) or iMarastonl d2005l), compared to 
the defau l t calib ration from iBruzual & Ch ariot (20QH) in 
iBell et all d2003l) makes no difference, since ellipticals are 
almost all relatively old and the model differences focus on 
young stellar populations. Assuming, therefore, that there 
is not some dramatic error in the stellar population models 
or the observed effective radii, then there are only a limited 
number of possible explanations for the FP "tilt." 

Non-Homology: Technically, this refers to systems being 
not perfectly self-similar. This must be broken in any model 
in which there is a dependence or "tilt" in the physical mass 
within some radius and the stellar mass M». In practice, 
the more practical meaning of "homology" is that systems 
have the same ratio of dynamical mass to true enclosed to- 
tal mass within their effective radii; i.e. dynamical mass is a 
good proxy for true mass without a systematic dependence 
on other galaxy properties. In other words, M true = kR e a 2 /G 
where k w constant (M tme oc R e a 2 /G is expected on dimen- 
sional grounds; k is the integral argument which depends on 
e.g. the details of profile shape and kinematics). There are two 
general possibilities: 

(a) Kinematic Non-Homology: In this case, the meaning 
of a changes with mass, either because of different contri- 
butions of rotationally supported components, or because of 
varying isotropy. Essentially all studies have found that this 
effect does not contribute to the FP tilt (see e.g.lGerhard et alJ 
| 2001[lRiciputi et alJ2005HCappellari et al.ll2007b iNipoti et al.l 
2008, and references therein); specifically that there is no cor- 
relation between FP residuals and orbital anisotr opy or rota- 
tion (w ithin the elliptical population). Moreover, Cio tti et all 
(1996) demonstrated that reproducing the observed tilt with 
such trends would require models with internally inconsistent 
or unphysical orbit structure. 

(b) Structural Non-Homology: Alternatively, the shape of 
the light profile could change sufficiently, such that e.g. even 
if two systems were both spherical with isotropic velocity 
dispersion tensors, the integral term relating M tme and Md yn 
is significantly different. (Note that this is not really inde- 
pendent of (a); the dominant effect of a change in the mass 
profile shape is generally to change the central velocity dis- 
persion, and correspondingly Md yn , rather than dramatically 
changing e.g. the physical meaning of R e , and differences in 
mass profile shape will generically require some kinematic 



Dissipation and the FP 



25 



non-homology.) There has been debate r egarding this possi- 
bility , as some observatio nal studies (e.g. iPrugniel & Simienl 
119971; iTruiillo et a"T]|2004l) have argued that the observed de- 
pendence of galaxy Sersic index on mass is sufficient to drive 
the tilt in the FP via structural non-homology. Others have 
noted, however, that integral models which allow for multi- 
ple components or do not explicitly assume a constant M/L 
with radius do not find s uch non-homology (e . g.lCarollo et aU 
[1991 iRixetail [19971; iGerhard et all l200~lt iGerhardl 120031; 
IPadmanabhan et al.ll2004t ICappellari et al.ll2006l) . Moreover, 
invoking structural non-homology can generally explain only 
differences in central velocity dispersion (the predicted FP 
tilt from such dynamical models disappears rapidly as a is 
measured at somewhat larger fractions of R e ); the fact that 
the FP tilt is not dramatically reduced using velocity disper- 
sions measured in larger aper tures (even out to R e ) rules out 
most classes of such models (iBender et al.lfl 994; Ciot ti et all 
1996t ISimien & Prugniellll997L 120021: lEmsellem et all 120041; 
Cap pellari et al Il2006t iNipoti et al .1 120081) . 

Our analysis is able to explain both sides of this debate. 
Fitting galaxy light profiles in a simplified manner, to a single 
Sersic index, does appear to yield a significant dependence of 
Sersic index on stellar mass, luminosity, or size. In Paper I, 
we demonstrate that this is true of our simulations - i.e. when 
analyzed in the same manner as the observations, they ap- 
pear to yield an identical dependenc e of Sersic in dex on mass 
to that found in IPrugniel & Simienl (119971) and ITruiillo et al.l 
(2004) (less massive systems, chosen in the same manner as 
those here to have higher typical dissipational fractions, have 
lower fitted Sersic indices over the same observational dy- 
namic range; see also Paper III) - even though (at fixed degree 
of dissipation), the systems are effectively homologous. Why, 
then, does this occur? 

The answer owes to the fact that a single Sersic index is not 
a physically robust description of the light profile of an ellip- 
tical over its entire extent, and it also fails to include the dark 
matter distribution. Such a fit mixes the outer dissipationless 
component of the remnant with the inner, compact starburst 
remnant. All the observed samples which appear to find a de- 
pendence of Sersic index on mass, when re-analyzed based 
on our physically motivated two-component decomposition, 
in fact reveal that the outer, dissipationless Sersic-like compo- 
nent is self-similar (see Paper II and Paper III); what changes 
with mass or radius is the mass fraction in the central dissi- 
pational component. Furthermore, because of the true multi- 
component nature of galaxy light profiles, fitting a single Ser- 
sic index yields different answers as a function of the dynamic 
range employed: for fixed observing conditions, this yields 
an apparent correlation of Sersic index with galaxy mass or 
size ( see e.g. IPadmanabh an et al1l2004t [Boylan-Kolchi n et all 
2005). Because of the role of e.g. the mass fraction at large 
radii and subtle differences in the second derivatives of pro- 
file shape near R e in driving a fitted Sersic index, the ob- 
served dependence of Sersic index on mass reflects a complex 
mix of real differences in galaxy merger history (the fact that 
low-mass systems have undergone more dissipational merg- 
ers, and that high-mass systems are increasingly likely to have 
undergone subsequent re-mergers, which will conserve the FP 
but scatter some small stellar mass to larger radii, leading to 
larger fitted Sersic indices) and observational differences in 
dynamic range, but little significant non-homology. 

We have demonstrated here (§ ISJ and in Paper II that dif- 
ferences in the degree of dissipation (i.e. in the strength of 
the central "extra light" component) - while capable of driv- 



ing substantial differences in the, fitted Sersic index when the 
galaxy is fit to a single Sersic law - do not in fact drive tra- 
ditional non-homology. They are not a large enough fraction 
of the total mass to dramatically alter the structural integrals. 
This should not be surprising: the dependence of M^/ M& yn 
on Sersic index (in models where the mass distribution fol- 
lows a pure Sersic law) is primarily driven by the behavior 
at large radii - large Sersic indices asymptotically approach 
the power-law behavior I(r) oc r" 2 , which implies a divergent 
mass at large r, and an implied effective radius R e — > oo and 
rapidly increasing central velocity dispersion. This is pre- 
cisely where observed profiles are in fact self-similar (in a 
mean sense; there is considerable variation object-to-object), 
with at most a weak dependence on formation history (i.e. 
median difference An s ~ 1 for those which have undergone 
major "dry" mergers; much less than is needed to explain the 
FP tilt), as we demonstrate in both observations and simula- 
tions in Paper II and Paper III. 

Further, at these radii, the mass density is increasingly dark 
matter-dominated: the structure integrals only depend signif- 
icantly on the Sersic index n s if the dark matter profile is also 
described by the same Sersic index - i.e. if the total M/L is a 
constant function of radius, whereas essentially all kinematic 
data at large radii indicates this is not the case (dark matter is 
required to explain the kinematics at large radii). As we have 
shown, and as demonstrated in observational kinematic mod- 
eling (e.g. IPadmanabhan et aU 120041; ICappellari et aT1l2006l) 
and lensing studies (Bolt on et al.l l2007, 2008), the dark mat- 
ter distribution is relatively insensitive to the details of the 
stellar mass distribution shape, and therefore the actual non- 
homology driven by even large changes in the stellar n s is 
much less than would be calculated assuming the dark matter 
followed the same profile (if it were not so, the observed large 
dispersion in n s values would necessarily yield much larger 
scatte r about the FP than is observed; see e.g. iBertin et all 
120021) . 

The implication that a change in n s should lead to signifi- 
cant structural non-homology is therefore misleading: the ob- 
served dependence of n s on mass owes mainly to issues of 
finite dynamic range and varying dissipational fractions, none 
of which give rise to significant structural non-homology. 
However, extrapolating to large radii based on these estimated 
n s values and the incorrect assumptions (1) that a single Sersic 
law is a physically meaningful parameterization of the galaxy 
light profile, and (2) that dark matter traces the stellar mat- 
ter in a strict one-to-one fashion, implies (incorrectly) some 
structural non-homology. 

Moreover, many studies indicate that dynamical mass is 
proportional to true enclosed mass (without requiring any 
homology assumption); i.e. that there cannot be a large con- 
tribution to tilt from any form of traditional non-homology, 
whether kinematic or structural. There ar e now indepen- 
dent l ines of supp ort for this conclusio n: ICappellari et al 
(120061) (see also Ivan der Marel 1199 It iKronawitter et al 



2000; Haring & Rixl 120041) estimate true enclosed masses 
within R e based on two and three-integral modeling, from 
two-dimension al velocity maps o f o bserved ellipticals. 
Alternatively, iBolton et aU (I2007L 120081) measure strong 
lensing gravitational masses. In both cases, the authors find 
M true oc Md yn , without any significant dependence on mass - 
in other words, the FP is unchanged regardless of whether 
the true total mass enclosed in R e is used, or whether the 
dynamical mass proxy Md yn is used. It appears that at most, 
at the ~ 3 a limits on the observed M true -Md y „ relation (and 
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based on the fitted M true -M* relations in Figure |2TTi. these 
traditional forms of non-homology may contribute ~ 1 /4 the 
observed tilt. 

Genuine Change in M/L: The observations therefore im- 
ply that the FP must reflect a genuine physical difference as 
a function of mass: namely, that low-mass ellipticals have a 
higher ratio of stellar or baryonic mass to total (baryonic plus 
dark matter) mass within the stellar R e , relative to high-mass 
ellipticals. All means of achieving this end can be classified 
into one of three categories: 

(a) Varying Global Baryon Fractions: One could imagine 
that the stellar and halo mass distributions of different ellipti- 
cals are separately self-similar, but that the total ratio of stellar 
to halo mass changes as a function of galaxy mass. In other 
words, there are no structural changes implied: one decreases 
the halo mass relative to galaxy mass in lower-mass systems. 

While this may seem plausible at the highest galaxy masses 
(it is well established that in bright clusters, the total stellar 
to dark matter mass ratio decreases with mass), this is rele- 
vant over much larger scales than the stellar R e of the central 
galaxy. Moreover, over most of the mass range of interest, 
the known trends in global stellar to dark matter fractions are 
opposite to that needed to explain the FP tilt. 

At masses below ~ M», where the FP is observed to be con- 
tinuous (recall, the observed FP and its tilt exte nd to systems 
with masses < 10~ 2 M*), all ACDM models (e.g.lConrov et alJ 
2006 ]IZentner et al.ll2005HZheng et al.ll2005t IVale & Ostrikerl 
2006t IShankar et alJ I2006t Ivan den Bosch et all 120071) and 
observational constraints (lEke et al.ll2004t lYang et aLl 12005b 
Mandelbau m et all 12006: Weinma nn et alJ 12006: Wan g et al.l 
2006) require that the ratio of total dark matter halo mass to 
stellar mass is higher in lower-mass systems - i.e. that star 
formation is less efficient in low-mass systems. This is con- 
trary to the effect desired here, and demonstrates that the FP 
tilt cannot owe to a simple global change in baryon fraction. 
Indirectly, however, the cosmological trend of stellar to dark 
matter mass is in fact important - this lower star formation ef- 
ficiency in low-mass systems means that they have larger gas 
fractions when they undergo mergers, which we show does 
give rise to the FP tilt. 

If the global baryonic mass ratio is unchanging at small 
radii (or changes in the opposite sense needed to tilt the FP), 
then the only other possibility is that the size/shape of the halo 
and stellar distributions vary relative to one another: i.e. one 
of the two components is made more or less compact, relative 
to the other, altering the ratio of stellar to total mass within the 
stellar R e . Recall, the radii of interest are small relative to the 
halo virial radii, and contain only a small fraction of the total 
halo dark matter mass, so regardless of the total halo to stellar 
mass ratio, changing the compactness of one component rela- 
tive to the other can make a significant difference. There two 
possibilities for this relative contraction/expansion: 

(b) Baryons are Fixed, Halos Contract: In this scenario, 
the galaxy stellar mass distributions are "scale free" (i.e. do 
not change owing to external factors), but the dark matter ha- 
los are less compact in low-mass ellipticals, so that within 
the stellar R e (i.e. the central regions of the halo) the total 
dark matter mass fraction enclosed is lo wer. This is quickly 
ruled out: all cosmological models (e.g. iNavarro et al.l 1996 ; 
Bullock et al1l200lHWechsler et al.l l2"002t iDolag etafll2004 ; 



find that halo concentration is a weak, decreasing function of 
galaxy mass, the opposite of the effect desired. 

Furthermore, if this were the case, one would expect to 
see it in progenitor disks, as well (the halos are insensitive 
to the morphological transformation of their central galaxies) 
- however, the baryonic Tully-Fisher relation and constraints 
on the dark matter halos of disks (e.g.|Persic & Saluccil 1990 ; 



on tne dark matter nalos or disks (e.g. Fersic &l ialucci 
iPersic et al.l fl996b: Bel l & de Jongl 120011: iMcG augh 
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Kuhlen et al. 2005; Neto et al. 2007) and observations from 



e.g. weak lensing and X-ray mass measurements dBuote et al.l 
120071: iSchmidt & AUenll2007T: iComerford & Nataraianll2007l) 



and references therein) reflect the expected cosmological 
trends (namely, baryon fractions and halo sizes scaling as pre- 
dicted in ACDM, in a weaker and opposite sense from that 
necessary if the FP tilt were to be explained in this manner). In 
short, effects (a) and (b) would predict that disks should fol- 
low the same FP scalings as ellipticals (modulo possible nor- 
malization offsets) - while in fact, they obey different scalings 
with, in many cases, an opposite qualitative sense (e.g. their 
size-mas s, velocity-mass, surface brightness-mass, and FP 
scalin gs; iFaber & Jacksonll 19761; iKormendvll 19891; IShen et alJ 
120031) . 

Figure [24] demonstrates this explicitly. We show the stellar 
size-mass relation of ellipticals ( from Figure \T% , compared 
to that of disks (from lShen et alJ (12003), but for our purposes 
here different sources agree well). Modulo a small normal- 
ization offset (owing to the profile shape), ellipticals would 
necessarily obey the same correlation if they were formed in 
purel y dissipationless mergers (see also ICiotti & van Albada 
2001). Obviously, ellipticals at low mass (where we empir- 
ically estimate and theoretically expect dissipation to be im- 
portant) are much more compact in their stellar distributions 
than disks. We also compare the FP of both types of ob- 
jects; i.e. dynamical mass Md yn (measured within R e of the 
stellar light) versus stellar mass. We estimate the relation for 
disks b ased on th e data and best-fit relations in Bell & d e Jon 
(120011) (see also IPersic & Saluccil 1 19901 IPersic et alJl!996 
who construct a similar correlation). Note that, for our pur- 
poses here, it makes little difference whether we plot the bary- 
onic or stellar mass (there is almost no difference for ellipti- 
cals, owing to their small gas content, and for disks, the re- 
lations fall within the quoted uncertainties in either case). It 
also makes no significant difference whether we use the same 
virial constant k to estimate Md yn for both disks and ellipticals, 
or attempt to make some correction for either profile shapes or 
the use of a velocity dispersion as opposed to a circular veloc- 
ity (the difference is small, a factor < 2, and we are ultimately 
interested in the qualitative scalings). 

It is clear (in agreement with previous work) that the scaling 
of Mdyn/M* in disks is opposite that of ellipticals: there is 
either no tilt or inverse tilt in the FP (i.e. M^jM^ is the same 
or higher in low-mass disks). For ellipticals where we have 
kinematic data as a function of radius, we can test whether the 
difference in the disk and elliptical scalings owes to effects 
(a) or (b); we do so by evaluating their enclosed mass Md yn 
not at the observed stellar effective radius of the spheroid, but 
at the radius of an equivalent disk (i.e. at the mean R e of a 
disk of the same stellar mass). The resulting trend of Md yn 
versus is indistinguishable from that of observed disks - 
i.e. by considering elliptical properties at their equivalent disk 
radii (radii they would have in the absence of dissipation), we 
effectively remove the tilt of the FP, and recover the observed 
correlations of disks. 

At the radius of a disk of similar mass (equivalently, at the 
radius the elliptical would have, if it were the product of a 
purely dissipationless merger of stellar disks), ellipticals have 
the same enclosed total mass as equivalent disks. In other 
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F IG. 24. — Top: Mean stellar size-mass relation o f ellipt icals, from Fig- 
ure [19] compared to that of disks (from IShen et all f2003lV ). Middle: FP 
of spheroids and disks - i.e. ratio of Mj VI 1 (ev aluated at R e ) to M*. We 
show the observed ellipticals from Figure 1151 and the best-fit power law 
Mdyn oc M\ +a to both relati ons (with correspondin g uncertainties). Disk cor- 
relation is from the data in Bell & de Jong (2001). At a given stellar mass, 
(low mass) spheroids have more compact stellar mass distributions, and less 
total mass (Mi yn ) enclosed within their stellar S c - Bottom: Same, but measur- 
ing M t | yn for the same observed ellipticals (from middle) at the mean expected 
radius of an equivalent (similar stellar mass) disk (from top). The best-fit re- 
lation to these data is shown; it is indistinguishable from the relation for disks. 
At the radii of their equivalent disks, ellipticals and disks have the same en- 
closed total (dark matter plus stellar) masses. The distribution of dark matter 
is not significantly different - the smaller sizes and M^ yn of ellipticals must 
reflect a contraction of the baryonic material relative to the dark matter, as 
predicted to occur in dissipative mergers. 

words, at the same (equivalent dissipationless) radius, disks 
and ellipticals of the same mass have the same dark matter 
mass content and distribution. In general, it is observationally 
well-established that the FP correlations of ellipticals become 
less distinct from those of disks as their properties are mea- 
sured at larger radii. These correlations clearly rule out sce- 
narios (a) and (b) above: if case (b) were true, disks and ellip- 
ticals should obey a similar stellar size-mass relation, and at 
the R e of an equivalent disk, ellipticals should still have much 
lower Md yn than spirals. If case (a) were true (i.e. both stel- 
lar and dark matter were more compact in ellipticals, but with 
lower total dark matter to stellar mass ratios), then we would 
again expect much lower Md yn at the equivalent disk radius in 
ellipticals. 



(c) Halos are Fixed, Baryons Contract: The only remain- 
ing possibility is our conclusion in this paper: namely, that the 
FP tilt arises because lower-mass spheroids have more com- 
pact stellar mass distributions, relative to their halos (equiv- 
alently, relative to what their stellar mass distribution would 
be in the absence of dissipation). We have demonstrated that 
this outcome is a natural consequence of dissipation in merg- 
ers - regardless of the initial scalings obeyed by progenitor 
disks, the fact that low mass disks are more gas-rich implies 
that, on average, there will be more dissipation in their merg- 
ers, yielding more compact baryonic remnants (while having 
little effect on the halo), and therefore increasing the ratio 
of stellar to dark matter mass inside the stellar R e in lower- 
mass systems. We further demonstrate that this reproduces 
precisely the observed tilt and scalings of elliptical properties 
with mass, while being consistent with all other observational 
constraints on the FP, including the "homology constraint," 
thatM tlue ocM dyn . 

We have also shown that dissipation is necessary in achiev- 
ing this. It has been known for some time that dissipation- 
less mergers cannot alter the phase-space density of ellipti- 
cals relative to their disk progenitors; consequently, ellipti- 
cals produced through dissipationless mergers obey the same 
scaling relations as their spiral progenitors (modulo normal- 
ization offsets), in stark contrast to the observed FP scalings. 
We demonstrate this explicitly: systems with the same dissi- 
pational fraction have the same ratio of total to stellar mass 
within R e ; i.e. do not internally exhibit any FP "tilt." Fur- 
thermore, we show that dissipation is the dominant factor de- 
termining the effective radii of ellipticals at fixed mass, even 
allowing for differences in formation and merger history - 
therefore if the explanation for the FP invokes any systematic 
change in elliptical sizes, it must involve dissipation. Only 
when the observed dependence of dissipation on mass is in- 
cluded is the observed tilt recovered. 

11.2. Additional Predictions 

We have extensively considered the role of dissipation in 
setting the FP tilt, effectively changing the ratio of Md yn /M„. 
To the extent that other properties also trace the degree of dis- 
sipation, we predict that these should similarly correlate with 
Md yn /M*. In Paper I and Paper II we develop an extensive set 
of predictions for elliptical properties that relate to the degree 
of dissipation in the spheroid-forming merger, and show how 
these relate to, e.g., the observed extra light (i.e. the tracer 
of the degree of dissipation in the spheroid forming merger- 
induced starburst). We refer to those papers for details and a 
large number of observational proxies of the degree of dissi- 
pation which can be used to further test the ideas herein. 

To the extent that the degree of dissipation in the original 
spheroid-forming merger reflects the gas fractions of the pro- 
genitor disks, it must also reflect the progenitor star forma- 
tion history. Broadly speaking, disks with more extended star 
formation histories would be expected to have larger gas frac- 
tions at the time of their merger, with younger stellar popula- 
tion ages and lower a-enhancement. If we ignore the effect of 
the merger on these stellar populations (a reasonable assump- 
tion if the system is observed at times significantly later than 
the merger, and if the mass fraction formed in the merger- 
induced starburst is not large, which are true for most of the 
observed systems of interest here), then these should be re- 
flected in the stellar populations of the ell iptical remnant. Pre- 
liminary observational comparisons from iGraves et al.l f2008) 
appear to support these predictions, and can provide powerful 
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independent tests and constraints for models of dissipation in 
spheroid formation: we therefore outline some relevant quan- 
titative predictions. 

Consider the following highly simplified toy model: iden- 
tical progenitor disks with initial gas fraction / gas = 1 fol- 
low an exponential, closed-box star formation history with 
time scale t, i.e. oc exp(-f/r), and merge at time t m , 
when the remaining gas is rapidly consumed in a central star- 
burst. The gas fraction at the time of the merger (and cor- 
respondingly, the dissipational fraction in the starburst) will 
be / sb = exp(-f m /r) (giving r = f,„/hi(l// sb )), and the mass- 
weighted mean formation time of the stars will be ff 0rm = 
t [1 -exp(-f m /r)] = f m (l -/ sb )/ln(l// sb ). For systems with 
a similar redshift range of their last major merger (similar t m , 
expected for systems of similar mass), their stellar formation 
times should therefore correlate with the starburst fraction / sb 
(both depending implicitly on the pre-merger star formation 
timescale r). 

We have shown in §[6] how M^n/M* is predicted to scale 
with the dissipational fraction / sb (this can be roughly ap- 
proximated as Mdyn/M* ~ (/ sb /0.2) -1 / 2 over the range of in- 
terest); combining the two, this yields an expected correla- 
tion between formation time and x = Mi yn /M* of the form 
f form = 0.5f„,(l -0.2x" 2 )/ln(2.24x). Observed at z = (i.e. 
with age ?h -fform), ellipticals with larger x = M& yil /M* at fixed 
mass (similar t m ) should be older - for a typical t m ~ lOGyr, 
systems with Md yn /M* w 2 are predicted to be ~ 2Gyr older 
than systems withMd yn /M* w 1. 

A shorter star formation history also implies higher a- 
enrichment in the progenitors. If we adopt the correlatio n 
for simple star formation models in Tho mas et al.l (2005), 
[a/Fe] w 1/5- 1/6 log A t (where Af is the star formation 
timescale for a Gaussian burst, but for our purposes here can 
be replaced by r modulo a conversion constant), and the scal- 
ings above, then we obtain the result (independent of the 
merger time t m ) that ellipticals with higher M^ yR jM^ should 
be more a-enriched. Specifically, objects with Md yn /M* « 2 
are predicted to have [a/Fe] values w 0.04-0.05 higher than 
systems withMd yn /M» w 1. 

Predictions for the absolute metallicities are more ambigu- 
ous (but see Paper II). In general, the trend of total/mean 
metallicity with mass will be dominated by the metallic- 
ity of pre-merger disks, which observations show (excluding 
the most gas-rich disks, where self-enrichment in the merger 
will dominate the final total metallic ity) trace a sim ilar mass- 
metallicity correlation to ellipticals (Gallazzi et al. 20051). In 
the absence of outflows or recycling, the metallicity would be 
the same for any systems with the same total accreted gas con- 
tent and stellar mass (with no dependence on the "starburst" 
content at fixed mass). However, to the extent that outflows 
in low-mass systems are responsible for the mass-metallicity 
relation (as is generally believed), then the detailed interplay 
between these outflows and merger-induced starbursts will be 
important. In broad terms, if outflow strengths and velocities 
are similar, then we expect dissipational star formation at the 
center of the galaxy to retain a higher metal content in com- 
parison to the same star formation in a more extended disk 
(since escape velocities from the galactic center and densi- 
ties leading to radiative losses in the outflows are higher). In 
this case, at fixed mass, more dissipational (lower Md yn /M», 
higher surface brightness) systems should have slightly higher 
metallicities than their less dissipational counterparts (and, as 
demonstrated in Paper I, this should be correlated with the ef- 



fects above - at fixed mass, the dependence of both quantities 
on dissipation should give rise to an inverse correlation be- 
tween metallicity and stellar population age or a-enrichment). 
Experimenting with e.g. different degrees of dissipation, out- 
flow strengths, and initial disk metalliticies in our simula- 
tions, we estimate this to be a relatively small (~ 0.1 dex) 
effect - not sufficient to dramatically effect the global mass- 
metallicity relation, but potentially visible in detailed studies. 

There is one important caveat here: systems of the same 
mass might also have had more gas-rich progenitors because 
they formed from mergers at very early times (i.e. had a differ- 
ent merger time t m , in the toy model illustrated above), making 
them older a nd more a-enriched. However, cosmological es- 
timates (e.g. Hopkins et al. 2008b) suggest this process is not 
dominant at a given stellar mass - i.e. systems of comparable 
stellar mass (and correspondingly similar total halo masses) 
tend to have similar merger histories. Specifically, the rela- 
tively large scatter in star formation history and disk gas frac- 
tions at fixed mass and redshift (a factor ~ 2 in / gas ) is larger 
than the scatter introduced by the combination of a scatter in 
formation times and the systematic cosmological evolution of 
disk gas fractions with time. Furthermore, systems with such 
early mergers will usually have multiple subsequent mergers 
at later times, so they will grow significantly in mass and have 
their effective radii substantially modified by these additional 
processes (such that they should and will be compared to dif- 
ferent systems at z = 0). 

Considering higher order effects, we demonstrate in Paper 
II that the strength of stellar population gradients is correlated 
with the degree of dissipation in the original spheroid-forming 
gas-rich merger, and show in Paper III that this holds even 
for remnants of subsequent gas-poor "dry" re-mergers. The 
most useful gradients in this sense are metallicity gradients 

- stellar age gradients and (especially) color gradients evolve 
strongly with time even in a fixed, passively evolving ellipti- 
cal (owing to the change in relative M/L for young and old 
stellar populations) and as such are more ambiguous, and gra- 
dients in a-enhancement are more sensitive to the gradients 
and overall star formation histories in the pre-merger disks. 
Metallicity gradients are, on the other hand, generally domi- 
nated by the degree of dissipation and imprinted in the gas- 
rich merger, and are not sensitive to the star formation history 
of the pre-merger disks, making them a more robust diagnos- 
tic for our purposes. At fixed mass, stronger gradients indicate 
more dissipation, and so we predict that, at fixed mass, ellip- 
ticals with higher Md yn /M» should have weaker metallicity 
gradients (see Paper II). 

11.3. Summary 

We have demonstrated from observations that the tilt of the 
FP owes to differential degrees of dissipation as a function of 
mass. Lower mass disks are more gas-rich, so their mergers 
are more dissipational: a larger fraction of the remnant mass 
is formed in a dissipational, merger-induced compact central 
starburst in the final stages of a major merger. This yields 
a remnant with a more compact stellar mass distribution, i.e. 
smaller R e relative to their progenitor disks, in lower-mass el- 
lipticals. The dark matter distribution is only weakly affected 

- implying that the stellar distribution in low-mass ellipticals 
is more compact, relative to the dark matter, than in equiva- 
lent disks or higher-mass ellipticals. Consequently, relatively 
less dark matter mass is enclosed within the stellar effective 
radius in low-mass ellipticals, so the ratio of enclosed mass 
Md yn /M* is an increasing function of mass (Md yn /M» oc M"). 
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This is the "tilt" in the FP. Given the observed, quantitative 
dependence of gas fractions on mass, the tilt is predicted to be 
exactly that observed, a w 0.2. 

Using a new empirical method, we robustly estimate the 
amount of dissipation involved in the formation a given ellip- 
tical. Specifically, with data of sufficient quality, we separate 
the observed surface brightness profile into an outer, violently 
relaxed component, which was established in a dissipation- 
less manner, and an inner "starburst remnant" or "extra light" 
component. We demonstrate in Paper I, Paper II and Paper III 
that observations of both evolved ellipticals and recent merger 
remnants support the proposal that this compact nuclear mass 
component a good proxy for the true mass formed in a dis- 
sipational starburst. Using this proxy, we show that the ob- 
served sizes of ellipticals, at fixed mass, depend strongly on 
the degree of dissipation involved in their formation (more 
so than even e.g. the number of mergers in their formation 
history). Correspondingly, we show that the ratio of total 
to stellar mass within the stellar effective radius, M&yn/M*, 
is a function of dissipation, both globally and at fixed mass 
(in the sense that elliptical sizes and M&yn/M* are decreasing 
functions of the amount of dissipation). These observed de- 
pendences are highly significant (P nu n < 10~ 7 ). Motivated by 
this, we show that by removing the mean systematic depen- 
dence of dissipation on mass, we can empirically remove the 
tilt of the FP. Considering ellipticals with the same dissipa- 
tional extra light fractions, we show that they obey a relation 
Mdyn oc M» (i.e. ellipticals with the same extra light content 
have the same ratio of dynamical to stellar mass within R e , 
independent of mass). 

In the proposed di s sipational models of the FP ( Bekk ij 19981; 
Onorb e et al.1 120051: IRobertson etafl l2006ri iDekel & Coxl 
2006), the tilt of the FP, and its projected correlations (e.g. the 
steepness of the stellar size-mass correlation of ellipticals rel- 
ative to that of disks), arise because low-mass disks are more 
gas-rich, and therefore low-mass mergers and ellipticals will 
have (on average) systematically higher degrees of dissipa- 
tion, and therefore smaller (relative) R e and M^/M*. If we 
consider e.g. simulations that obey the observed correlation 
between disk gas content and mass (as opposed to being dis- 
sipationless, or having all the same gas fractions independent 
of mass - neither of which is consistent with observations), 
then the FP predicted has exactly the observed tilt. Equiva- 
lently, the observed mean dissipational fractions of ellipticals, 
as a function of mass, agree well with the observed gas frac- 
tions of progenitor disks of the same masses, over the redshift 
range z ~ — 2. In other words, dissipation is both necessary 
and sufficient to explain the FP tilt and differences between 
disk and elliptical scaling relations. 

To our knowledge, this is the first explicit observational 
test of these theoretical models. We further demonstrate that 
other mechanisms cannot be responsible for the majority of 
the FP tilt. For example, observations have demonstrated 
that the "homology assumption," namely that M tme oc a 2 R e , 
is valid, and we show that simulations predict this - the ho- 
mology breaking introduced by dissipation is negligible. In 
other words, the FP tilt reflects the ratio of stellar to true mass 
enclosed within R e : M tme /M* oc M" (this ratio is not con- 
stant), rather than an "apparent" effect. We also show that, 
if we measure elliptical properties at the radius of an equiv- 
alent disk, the tilt of the FP is removed: within the radii of 
disks of the same mass, ellipticals have the same total and dy- 
namical masses. That is, the dependence of Md yn /M* cannot 
be driven primarily by changes in the dark matter distribution 



at fixed baryonic properties, nor by changes in the total dark 
matter to stellar mass ratio (integrated over the entire halo). 
The variation in Md yn /M* within R e must predominantly re- 
flect the change in size of the baryonic component: low-mass 
ellipticals have much more compact stellar distributions than 
similar-mass disks, and therefore have less enclosed dark mat- 
ter within that stellar R e , as predicted in dissipational theories. 

Together, these observational tests represent an important 
vindication of the "merger hypothesis," that ellipticals are 
formed by the gas-rich mergers of disk galaxies, and mod- 
els for the origin of the FP in dissipational m ajor merg- 
ers (IRobertson et all l2006bt IDekel & Coxl l2006j) . We ex- 
plicitly demonstrate that, regardless of subsequent gas-poor 
(spheroid-spheroid or "dry") re-mergers, the location of sys- 
tems with respect to the FP and e.g. elliptical size-mass and 
velocity dispersion-mass relations is primarily determined by 
the amount of dissipation involved in their formation: i.e. the 
gas content involved in the original, spheroid-forming merger. 
Gas rich mergers cannot be ignored in the formation of ellipti- 
cals. Not only have we demonstrated that the FP is consistent 
with the merger hypothesis, but that (given the systematic de- 
pendence of disk gas fractions on mass), a FP tilted in the 
manner observed is a necessary prediction of the theory. 

We have also shown that elliptical sizes, inferred dissipa- 
tional fractions, and the FP are completely consistent with 
the formation of ellipticals in mergers of disks with similar 
properties (sizes, gas fractions, dark matter halo masses and 
sizes) to those observed in low- redshift (z ~ 0- 1) disks. In 
other words, the sizes of ellipticals and their FP correlations 
do not require elliptical progenitors to be more compact than 
observed, low-redshift disks. Dissipation is sufficient to ex- 
plain the differences in their densities and sizes. The fact 
that, within the radius of an equivalent disk, ellipticals obey 
the same correlation between total and stellar mass actually 
implies that their dark matter halos (and presumably those 
of their progenitors) are not significantly more compact than 
those of low-redshift disks. 

This is important for the viability of the merger hypothesis, 
given the observations indicating t hat disk and halo sizes do 
not evolve strongly with redshift (Ravindranath et al j l2004t 
iBarden et all 120051: iTruiillo et all l2006t IZirm et all 12007b . 
This is not to say that elliptical sizes might not evolve with 
redshift (which is easil y possible if e.g. disk gas fractions sys- 
temati cally evolve; see lKhoch far & Silk l2006tlHopkins et al.l 
2007a), nor that ellipticals all formed at low redshift (indeed, 
if the disk size evolution is weak, then ellipticals form could 
rapidly at relatively early times and still resemble the products 
of low-redshift disks). However, it does imply that exotic pro- 
genitors - progenitors not found in the local universe - are not 
required to explain the observed correlations, surface bright- 
ness profiles (see Paper II and Paper III), or kinematics (see 
ICox et al]|2006bh of typical local ellipticals. 
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